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Abstract 

We construct a new class of positive solutions for the classical elliptic problem 
Au-u + uP = 0,p>2, in 
We establish a deep relation between them and the following Toda system 
c^/j' = (Jj-i-fj - ef^-fj+i in R, j = 1, . . . , /c. 

We show that these solutions have the approximate form u{x, z) ^ ~ fji^)) 

where w is the unique even, positive, asymptotically vanishing solution of ui" — w + 
Kj^ = in R. Functions fj{z), representing the multiple ends of 7i(x,z), solve the 
aforementioned Toda system, they are even, asymptotically linear, with 

/o = -oo < /i <C • • • <C A- < fk+i = +00. 

The solutions of the elliptic problem we construct have their counterpart in the 
theory of constant mean curvature surfaces. An analogy can also be made between 
their construction and the gluing of constant scalar curvature Fowler singular metrics 
in the sphere. 
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1 Introduction and statement of main results 

This paper deals with the classical scmilinear elliptic problem 

Au~u + uP = Q, u>0, inM^ (1.1) 

where p > 1. Equation (jl.ip arises for instance as the standing- wave problem for the 
standard nonlinear Schrodingcr equation 

typically p = 3, corresponding to that of solutions of the form ip{y,t) = u(y)e~**. It 
also arises in nonlinear models in Turing's theory biological theory of pattern formation 
[39] such as the Gray-Scott or Gierer-Meinhardt systems, [16l[T5j. The solutions of 
(jl.ip which decay to zero at infinity are well understood. Problem (|l.ip has a radially 
symmetric solution WN{y) which approaches at infinity provided that 



1 <p < 



if AT > 3 
+00 a N ^ 1,2, 



see [3H1 13] ■ This solution is unique [H] , and actually any positive solution to p.ip which 
vanishes at infinity must be radially symmetric around some point [Tl] . 

Problem (|l.ip and its variations have been broadly treated in the PDE literature 
in the last two decades. These variations are mostly of one of the two types: (jl.ip 
is changed to a non-autonomous problem with a potential depending on the space 
variable; or p.ip is considered in a bounded domain under suitable boundary conditions. 
Typically, in both versions a small parameter is introduced rendering (|l.ip a singular 
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perturbation problem. We refer the reader to the works [Tll^lSligifTUlfnifnifTTl 
[TSl [23l [25l [26l [34l [32l [33] and references therein. Many constructions in the hterature 
refer to "multi-bump solutions" , built by perturbation of a sum of copies of the basic 
radial bump wjy suitably scaled, with centers adjusted in equilibrium under appropriate 
constraints on the potential or the geometry of the underlying domain. 

Much less is known about solutions to this equation in entire space which do not 
vanish at infinity (while they are all known to be bounded, see |36|). For example, 
the solution wn of in induces a solution in K.'^+^ which only depends on TV 
variables. This solution vanishes asymptotically in all but one variable. For simplicity, 
we restrict ourselves to the case N = 2, and consider positive solutions u(x, z) to 
problem p.ip which are even in z and vanish as |a;| +oo, namely 

u{x, z) = u{x, —z) for aU (x, 0) e (1.2) 

and 

lim u{x, 2) = for all 2 e M. (1.3) 

|a:| — »+oo 

A canonical example is thus built from the one-dimensional bump wi, which we 
denote in the sequel just by w, namely the unique solution of the ODE 

w"-w + wP = 0, w > 0, inM, (1.4) 
w\Q) = 0, w{x) ~> as |x| +00, (1.5) 

corresponding in phase plane to a homoclinic orbit for the equilibrium 0. Using this 
function we can define a family of solutions u of equation p.ip with the properties (|1.2p - 
(jl.3[) setting u(x, z) := 'w{x — a)^ a e R. By analogy with the above terminology, we may 
call these solutions " single bump- lines" . A natural question is whether a solution that 
satisfies (|1.2p - (|1.3|) and which is in addition even in x must equal w{x). The solution w 
of was found to be isolated by Busca and Felmer in 0] in a uniform topology which 
avoids oscillations at infinity. On the other hand, a second class of solutions which are 
even both in z and x was discovered by Dancer in [7j via local bifurcation arguments. 
They constitute a one-parameter family of solutions which are periodic in the z variable 
and originate from w{x). Let us briefly review their construction: we consider problem 
(jl.ip with T-periodic conditions in z, 

u{x,z + T)^u{x,z) for all (x, z) e (1.6) 

and regard T > as a bifurcation parameter. The linearized operator around the single 
bump line is 

L{(j)) = 0^, + (f>^^ + (pwP-^ - 1)<I). 
It is well known that the eigenvalue problem 

+ (pw'^"' - 1)0 = A0, (1.7) 

has a unique positive eigenvalue Ai with Z{x) a positive eigenfunction. We observe that 
the operator L has a bounded element of its kernel given by 



Z{x) cos(-\/Ai-z), 
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which turns out to be the only one which is even, both in x and z variable, and in 
addition T ~ ^S=-periodic in z. Crandall-Rabinowitz bifurcation theorem can then be 

V -^1 

adapted to yield existence of a continuum of solutions bifurcating at this value of T, 
periodic in z with period Ts = + 0{S). They are uniformly close to w{x) and their 
asymptotic formula is: 

ws{x,z) = w{x) + 6Z{x) cos{^/%z) + 0((5^)e"l''l . 

We refer to the functions wg in what follows as Dancer solutions. 

The purpose of this paper is to construct a new type of solutions of p.ip in that 
have multiple ends in the form of multiple bump-lines, and satisfy in addition properties 
(|1.2[) - (|1.3[) . To explain let us consider the function 

k 

(x) = w{x — Qj). 

3=0 

While for numbers ai ^ a2 ^ ■ ■ • <C , lu* appears to be a very good approximation 
to a solution of the ODE (|1.4p . the only solutions of (|1.4p which go to zero at infinity 
must be single translates of w(x). Our main result in this paper is that, in the two 
dimensions, there exist solutions u{x, z) with a profile of this type for each fixed z. Of 
course in this case the numbers aj must be replaced by non-constant functions fj{z). 
In addition, it will be necessary to use as basic cells not just the standard bump-lines, 
but rather the wiggling bump lines found by Dancer. 

Thus, what we actually look for is a solution u{x, z) which is close, up to lower order 
terms, to a multi bump-line of the form 

fc 

w*{x,z) = ^ws^{x - fj{z),z), (1.8) 
for suitable small numbers 5j and even functions 

h{z) < /2(z) < ••• < fk{z), 

which have uniformly small derivatives. The functions fj cannot be arbitrary and they 
turn out to satisfy (asymptotically) a second order system of differential equations, the 
Toda system^ given by 

c^/j' = ef^-^-f^ - ef^-f^+^ in M, j = 1, . . . , fc, (1.9) 

with the conventions /o = — oo, Jk+i = +oo, where Cp is an explicit positive constant 
that will be specified later. In agreement with the symmetry requirement (|1.2p we 
consider even solutions of system (|1.9|) . namely 

fj{z)^fj{-z) forallzeR. (1.10) 

We observe that for an even solution f = (/i, . . . , fk) of this system, function defined 

by fell 

ia = ifal,---,fak), faj{z) := fij{az) + {j ^ ) log - (1.11) 
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is also an even solution of the system. Our main result in this paper asserts that given 
such an f whose values at z = are ordered, for all sufhciently small a there exists a 
multi bump-line solution of approximate form (jl.8|) with f replaced with f^j and suitably 
chosen small numbers Sj dependent on a. Thus, we fix numbers 



fe 



ai < 02 < ■ • • < flfe, ^ai=0, (1-12) 



i=l 



(as we will see shortly the latter condition can be assumed without loss of generality) 
and consider the unique solution f of system (jl.9p for which 

/j(0) = aj, /j(0) = 0, j = l,...,/c, (1.13) 

and their associated scalings fa. As an explicit example, we immediately check that for 
k ~ 2 we have 

-l(log(2A-e,)-.og ' ), 

1 / 1 \ (1^"> 

AW = -(log(2A-c,)-l„g^_;^), 

where 



ga2 — oi 

As we will see later, in the case of a general k the functions fj are asymptotically linear: 
the limits Vj = f[j{+oo), exist and 

k 

vi<U2<--- <Vk, ^Uj^Q. (1.15) 



For instance, for fc = 2 we have V2 = —vi — 5 y e°2-''°i according to formula (|1.14|1 
Besides, we have (globally) for a small 

.fal{z) < /q2(z) < • • • < .fak{z), f'aj{ + <=^) = '^j": 



and 



U,{z) = v,a\z\+h, + (j-:^)logl + 0(e-^"N), 



for certain scalars hj and ?9 > 0. These are standard facts about the Toda system that 
can be found for instance in |20| . Thus, each of the multiple ends of u(x^ z) is a bump- 
line that is nearly straight but bent, with an angle slightly distinct than the angles of 
the other ends. The Toda system is a classical model describing scattering of k particles 
distributed on a straight line, which interact only with their closest neighbors with a 
forces given by a potential depending on the exponentials of their mutual distances. Here 
the z variable is interpreted as time. In this context, corresponds to a setting in which 
the particles starting from the rest, scatter at slightly different, nearly constant small 
velocities whose average is zero. The latter fact follows the identity X]j=i faji^) = 0, 
which also implies conservation of the center of mass used in (|1.12p . 



Our main result is: 



Multiple-end solutions of autonomous elliptic problems 



6 



Theorem 1.1 Assume that N = 2 and p > 2. Given k > 2 and numbers aj as in 
(ji. jgp . for any sufficiently small number a > 0, there exists a solution of equation 
which satisfies conditions <\1.2\ - ^1.S\ . and that has the form 

k 

Ua{x, z) = ^W5^{X - faj{z), z) (1 + o(l)). 

i=i 

Here is the scaling 1^1.11^ of f, the unique solution of satisfying I\1.1S^ . and 

o(l) ^0, Sj ^0 as 0. 

Remark 1.1 By no means we have intended to state Theorem II. II in its most general 
form. For instance, the same result holds for the more general problem 

Au + g{u) = Q, u>0, inR2. (1.16) 

where g is of class C^, ^(O) = 0, g'{0) < 0, g{s)ds = for some c > 0. In this case 
a homoclinic solution w analogous to that of (|1.4p - (|1.5p exists. On the other hand, for 
problem (|l.ip in M^, > 3, wc may take as the basis of the construction the unique 
radially symmetric decaying solution wn-i{x) of (jl.ip in M^"^, provided that p < 
for which Dancer solutions wn_i s{x) are equally available. In such a case we look for 
solutions close to a function of the form 

k 

^^WN-l,Sjixi,X2, ■ ■ ■ ,XN-2,XN-1 - fa]{z),z), 

which are radial in the first N ~ 2 variables. In both extensions the necessary changes 
in the constructions are straightforward, so that for simplicity we only consider here 
the case iV = 2 for pTTjl . 

Remark 1.2 The structure of equation (jl.ip makes it natural to look for solutions u 
which are even both in x and z variables. It turns out that if the numbers Oj in p.l2p 
satisfy the symmetry requirement 

Oj = -flfe+l-j, j = 1, . . . 

then the solutions fo, satisfy similarly 

faj ^ -fak+i-], for all j = 1, . . . , fc. (1.17) 

Indeed, in such a case the construction of the solutions in Theorem 11.11 vields that 
they are even in both variables: 

Uct(— .T, z) — Ua{x, z) = Ua{x, — z) , for all {x, z) G M^, (1-18) 

In particular for k even this solution satisfies 

lim Ua{x,z)~0, for all z e M, lim Ua{x,z)^Q, for all a; e M. (1-19) 

X — *±oo z — »±oo 

By the well known result by Gidas. Ni and Nirenberg [14j a positive solution of equation 
(|l.ip that satisfies (|1.18p and the limit conditions p.l9p uniformly must be radially 
symmetric around the origin. Theorem 1 1 . 1 1 shows that uniformity cannot be relaxed in 
this classical result. 



Multiple-end solutions of autonomous elliptic problems 



7 



One of the striking features of the existence result in Theorem 1 which is a purely 
PDE result, is that its counterparts can be found in geometric framework. Indeed, 
there arc many examples where correspondence between solutions of (jl.ip and those 
of some geometric problem can be drawn. To illustrate this, we will concentrate on 
what is perhaps the most spectacular one: the analogy between the theory of complete 
constant mean curvature surfaces in Euclidean 3-space and the study of entire solutions 
of (jl.ip . For simplicity we will restrict ourselves to constant mean curvature surfaces in 
R"^ which have embedded coplanar ends. In the following we will draw parallels between 
these geometric objects and families of solutions of 

Embedded constant mean curvature surfaces of revolution were found by Dclaunay 
in the mid 19th century They constitute a smooth one-parameter family of singly 
periodic surfaces Dr, for r £ (0, 1], which interpolate between the cylinder Di = S^{1) x 
R and the singular surface Dq '■= limT-*Q D^, which is the union of an infinitely many 
spheres of radius 1/2 centered at each of the points (0,0, n) as n e Z. The Dclaunay 
surface Dr can be parametrized by 

Xr{x, z) = {^{z) coax, ip{z) sinx, tp{z)) £ Dt <Z M.^ , 

for (x, z) £ R X IR/27rZ. Here the function ip is smooth solution of 



(^') 

and the function is defined by 



/\2 



As already mentioned, when r = 1, the Dclaunay surface is nothing but a right 
circular cylinder Di — S^{1) x M, with the unit circle as the cross section. This cylinder 
is clearly invariant under the continuous group of vertical translations, in the same way 
that the single bump-line solution of (|l.ip is invariant under a one parameter group of 
translations. It is then natural to agree on the correspondence between 



The cylinder 
Di = S"! X M 



The single bump-line 
(x, z) I — > w{x) 



Let us denote by the unique radially symmetric, decaying solution of (jl.ip . Inspec- 
tion of the other end of the Delaunay family, namely when the parameter r tends to 0, 
suggests the correspondence between 



The sphere 



The radially symmetric solution 



(x,z) 



To justify this correspondence, let us observe that on the one hand, as the parameter r 
tends to 0, the surfaces Dj- resemble a sequence of spheres of radius 1/2 arranged along 
the cca-axis which are connected together by small catenoidal necks. On the other hand 
an analogous solution of (jl.ip can be built as follows. Let 5*^ = R x (0, R) and consider 
a least energy (mountain pass) solution in H^{Sb.) for the the energy 



\Vu[ 



Sr 



Sr 



Sr 
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for large i? > 0, which we may assume to be even in x and with maximum located 
at the origin. For R very large, this solution, which satisfies zero Neumann boundary 
conditions, resembles half of the unique radial, decaying solution 1^2 of Extension 
by successive even reflections in z variable yields a solution to p.ip which resembles a 
periodic array of radially symmetric solutions of (jl.ip . with a very large period, along 
the z-axis. While this is not known, these solutions may be understood as a limit of 
the branch solutions constructed by Dancer. 

More generally, there is a natural correspondence between 



Delaunay surfaces 

Dr 



Dancer solutions 

(x, z) I — > ws{x, z) 



To give further credit to this correspondence, let us recall that the Jacobi operator about 
the cylinder Di corresponds to the linearized mean curvature operator when nearby 
surfaces are considered as normal graphs over Di. In the above parameterization, the 
Jacobi operator reads Ji = ^ [d^ + 9^ + 1). In this geometric context, it plays the 
role of the linear operator defined in section 2 which is the linearization of (jl.ip about 
the single bump-line solution w. Hence we have the correspondence 



The Jacobi operator 



The linearized operator 



In our construction, the polynomially bounded kernel of the linearized operator L plays 
a crucial role. Similarly, the polynomially bounded kernel of the Jacobi operator Ji has 
some geometric interpretation. Let us recall that we only consider surfaces whose ends 
are coplanar, the Jacobi fields associated to the action of rigid motions are then given 

by 

{x,z)i — >cosa; and (x,z)' — > z cosx, 

which correspond respectively to the action of translation and the action of the rotation 
of the axis of the Delaunay surface Di . Clearly, these Jacobi fields are the counterpart 
of the elements of the kernel of L which are given by 

(x, z) I — > dx'w{x) and {x,z)i — >zdxw{x), 

since the latter are also generated using the invariance of the problem with respect to 
the same kind of rigid motions. 

Two additional Jacobi fields associated to Ji are given by 

(x, z) I — *-cosz and (x,z)i — > sinz, 

which are associated to the existence of the family Dt as r is close to 1, as can be easily 
seen using a bifurcation analysis, in a similar way that the functions 

and 



(x, z) I — > Z{x) cos{y/Xiz) 



(x, z) I — > Z{x) sin(vAiz), 



are associated to the existence of Dancer solutions when the parameter 6 is close to 0. 
These two bifurcation results have their origin in the fact that we have the correspon- 
dence between 



The ground state 1 of 
92 + 1 



The first eigenfunction Z(yx) of 



1 



■ pw 



p-i 
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both of them associated to negative eigenvahies. The fact that the least eigenvalue of 
these operators is negative is precisely the reason why a bifurcation analysis can be 
performed and gives rise to the existence of Delaunay surfaces close to Di or Dancer's 
solutions close to the bump- line w. 

With these analogies in mind, we can now translate our main result into the constant 
mean curvature surface framework. The result of Theorem 11.11 corresponds to the 
connected sum of finitely many copies of the cylinder ^^(l) x M which have a common 
plane of symmetry. The connected sum construction is performed by inserting small 
catenoidal necks between two consecutive cylinders and this can be done in such a way 
that the ends of the resulting surface are coplanar. Such a result, in the context of 
constant mean curvature surfaces, follows at once from [53]. It is observed that, once 
the connected sum is performed the ends of the cylinder have to be slightly bent and 
moreover, the ends cannot be kept asymptotic to the ends of right cylinders but have 
to be asymptotic to Delaunay ends with parameters close to 1, in agreement with the 
result of Theorem 11.11 

However there is a major difference. The Toda system which governs the level sets 
has found no analogy in the constant mean curvature surfaces. This is mainly due to 
the strong interactions in the elliptic equations. 

Another (older) construction of complete noncompact constant mean curvature sur- 
faces was performed by N. Kapoulcas [H] (sec also [23) starting with finitely many 
halves of Delaunay surfaces with parameter t close to which are connected to a cen- 
tral sphere. The corresponding solutions of (jl.ip have recently been constructed by A. 
Malchiodi in [24] . It should be clear that many more examples of solutions of p.ip can 
be found using the above correspondence and we shall return to this in a forthcoming 
paper. 

It is well known that the story of complete constant mean curvature surfaces in 
R-^ parallels that of complete locally conformally flat metrics with constant, positive 
scalar curvature. Therefore, it is not surprising that there should be a correspondence 
between these objects in conformal geometry and solutions of (jl.ip . For example, De- 
launay surfaces and Dancer solutions should now be replaced by Fowler solutions which 
correspond to constant scalar curvature metrics on the cylinder R x S"^^ which are 
conformal to the product metric dz^ + ggn-i, when n > 3. These are given by 

where z i — > v{z) is a smooth positive solution of 

(„')2_„2+n^„;^^ ^_2^2 
^ ' n n 

When r = 1 and v = 1 the solution is a straight cylinder while as r tends to the 
metrics converge on compacts to the round metric on the unit sphere. The connected 
sum construction for such Fowler type metrics was performed by R. Mazzeo, D. Pollack 
and K. Uhlenbcrk [31] (where it is called the dipole construction). N. Kapouleas' 
construction mentioned above was initially performed by R. Schoen |37j (see also R. 
Mazzeo and F. Pacard [28]). 
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2 Linear theory 

In this section we will consider the basic linearized operator. The developments pre- 
sented here are crucial for our paper later on. By w we will denote the homoclinic 
solution to u" ~ u + uP ~ such that ■w'{0) = 0. Let 

w) = + r{w)cb, fiw) = pwp-^ - 1. 

We recall some well known facts about Lq. First notice that Lq{wx) = i.e. has one 
dimensional kernel. Second we observe that 

1 u;(P+l)/2 

Ai = t(p-1)(p + 3), Z = 

correspond, respectively, to principal eigenvalue and eigenfunction. Except for Ai and 
A2 = the rest of the spectrum of is negative. As a eonsequcnec of these facts we 
observe that problem 

Lo{^)-e4'^h, (2.1) 

is uniquely solvable whenever ^ 7^ ±v^, for h E L^(R). Actually, rather standard 
argument, using comparison principle and the fact that Lq is of the form 

LoW = + rz(x)0, \qix)\ < Ce-^l"l, 

can be used to show that whenever h is for instance a compactly supported function 
then the solution of (|2.ip is an exponentially decaying function. 
Let us consider operator 

L(j) = Lo{(j)) + 4>zz, 

defined in the whole plane {x,z) G K^. Equation L{(p) ~ 0, has 3 obvious bounded 
solutions 

Wx{x), Z{x) COS \/Xiz, Z{x) sin{\/Xiz). 
Our first result shows that converse is also true. 
Lemma 2.1 Let (p he a hounded solution of the prohlem 

L{(t)) =0 in R2. (2.2) 

Then (j){x,z) is a linear combination of the functions Wx{x), Z{x) cos{^/\iz), and 
Z{x) sva{\f\\z) . 

Proof. Let assume that (/> is a bounded function that satisfies 

<t'.. + + {vw"'^ - 1)0 = 0. (2.3) 

Let us consider the Fourier transform of (/)(a;, z) in the z variable, 0(a;,^) which is by 
definition the distribution defined as 

(^(a;,-),Ai>R = (0(2;r),A>M = / H^^OKOd^, 
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where /i(^) is any smooth rapidly decreasing function. Let us consider a smooth rapidly 
decreasing function of the two variables Then from equation (|2.3p we find 

/ {^(x, •), - + {Pw""'^ - iWmdx = 0. 
Jm. 

Let ip{x) and fi(^) be smooth and compactly supported functions such that 

{ - 0} n Supp (m) = 0- 
Then we can solve the equation 

- + ipuj'^' - 1)^- = f^iOvix), x&R, 

uniquely as a smooth, rapidly decreasing function tp{x, ^) such that ip{x, ^) = whenever 
^ ^ Supp (/i) . We conclude then that 

{4){x, ■),^j,)ts.(p{x) dx = 0, 

so that for all x € M, < •), /.j >= 0, whenever {VAi, ^ 0} Supp {^) = 0, in 
other words 

Supp{(j>{x, •)) c {VaT, -yXi",o}. 

By distribution theory we find then that (j){x, •) is a linear combination, with coeffi- 
cients depending on z, of derivatives up to a finite order of Dirac masses supported in 
{VA]", —^/Xi, 0}. Taking inverse Fourier transform, we get that 

(j){x, z) = po{z, x) +pi{z, x) cos{\/Xi_z) +p2{z, x) sin{\/Xiz), 

where pj are polynomials in z with coefficients depending on x. Since (f) is bounded 
these polynomials are of zero order, i.e. Pj{z, x) = Pj{x), and the bounded functions pj 
must satisfy the equations 

■^o(po) = 0, Lo(pi) - Aipi 0, io(P2) - Aip2 = 0, 

from where the desired result thus follows, q 
Let B{(f)) be an operator of the form 

B{4>) = bidxx(t> + h2dxz(t> + h:ida:4> + hid,^(t) + 65^, 

where the coefficients bi are small functions. In the sequel we will denote b = (61, . . . , 65) 
and assume that 

.5 

||b|| = ^ \\b,\\oo + IIV61II00 + ||V52||oo < 5o, (2.4) 

where the small number i5o will be subsequently fixed. The linear theory used in this 
paper is based on a priori estimates for the solutions of the following problems 

+ L(0) = h, in R2. (2.5) 
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The results of Lemma |2 . II imply that such estimates without imposing extra conditions 
on (f) may not exist. The form of the bounded solutions of = suggests the 

following orthogonality conditions: 



'{x, z)wx{x)diJ,{x) = = J (j){x, z)Z{x)d^{x), for all z e R, (2.6) 

where d^{x) is a fixed measure in M absolutely continuous with respect to the Lebesque 
measure. In the sequel we will in particular consider dfj,{x) — p{x) dx where p is a 
compactly supported cut-off function, however our next result applies for a general 
dfi{x) as well. With these restrictions imposed we have the following result concerning 
a priori estimates for this problem. 



Lemma 2.2 There exist constants So and C such that if the hound \2.4-\ holds and 
h G L°°{JS?), then any hounded solution </> of prohlem (|i^.5p - (|i^.6'P satisfies 

||</'||oo < C\\h\\oo. 

Proof. We will argue by contradiction. Assuming the opposite means that there are 
sequences (/)„, /i„ such that 



|l<^«|loo = l, |l/ln|loo->0, 

and 

Bni^n) + L{(Pn) = K, inR\ (2.7) 



(/)n{x, z)wx{x)dfj,{x) = = / (pnix, z)Z{x)dii{x), for all z G M. (2.8) 
Here 

Let us assume that (a;„, z„) G is such that 

|0n(a;„,z„)| 1. 

We claim that the sequence Xn is bounded. Indeed, if not, using the fact that Lcj) = 
A(/) — + O(e~'^l^l)0 and employing elliptic estimates we find that the sequence of 
functions 

^„(X, z) = (j)n{Xn + X, Z„ + z), 

converges, up to a subsequence, locally uniformly to a solution (j) of the equation 

A^-0 = O, inR^ 

whose absolute value attains its maximum at (0,0), This implies = 0, so that Xn is 
indeed bounded. Let now 

(pnix, z) = <j)n{x, Zn + z) . 

Then 0„ converges uniformly over compacts to a bounded, nontrivial solution of 

L(^) = in R2, 
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{x, z)wx{x)dfi{x) = = / (l>{x, z)Z{x)dfi{x), for all z e R. 
Js. 

Lemma |2. II then implies = 0, a contradiction and the proof is concluded, q 

Using Lemma l2.2l we can also find a priori estimates with norms involving exponen- 
tial weights. Let us consider the norm 

||0||.,a^ l|e'^l^l+"l^l(/)||oo. 

where numbers cr, a > are fixed and will be subsequently adjusted. In the case a = 
we have the following a priori estimates. 



Corollary 2.1 There are numbers C and do as in Lemma \2.2\ for which, if \\h\\a,o < 
+O0, a S [0, 1), then a bounded solution (p of i2.5\ )- ^276\) satisfies 

11011^,0 + ||V0||,,o<C||/i||,,o. (2.9) 



Proof. Again we concentrate on estimates for the problem (|2.5p - (|2.6p . We already 
know that 

||0||oo<C||/l|U,O. 

We set ^ = (j)\\h\\^),. Then we have 

{L + B){^) = h, where II /i|U,o < 1, 
and also ||^||oo < C. Let us fix a number i?o > such that for x > Rq we have 

l-a2 



which is always possible since w{x) = 0(e~^'^'). For an arbitrary number p > let us 
set 

4>{x, z) = p[cosh(z/2) + e*"^] + Me'"'', 
where M is to be chosen. Then we find that, reducing 6q if necessary, 

{L + B){4>) < fora;>i?o. 

Thus 

{L + B){(j)) <h, forx>i?o, 

if 

M(l-o-2) 

->\\h\\afi = l. 

If we also also assume 

Me— > ||<^||oo, 

we conclude from maximum principle that (}> < (f>. Letting p ^ we then get by fixing 
M, 

4><Me-'''', fora;>0, 
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hence 

(p< MWhW^fiC-"'', forx>0. 
In a similar way we obtain the lower bound 

(/) > -il/||;i||„,oe-'"^, fora;>0. 

Finally, the same argument for x < yields 

muo < c\\h\\,,o, 

while from local elliptic estimates we find 

||V0||,,o < C\\h\\,,o, 

and the proof is concluded, q 

When a > in the definition of the norm j| • ||o-,a then we have the following a priori 
estimates. 

Corollary 2.2 There are numbers C, (5o o-s in Lemma \2/A and ao > for which, if 
\\h\\a.a < +00, CT € (0,1), a G [0, ao), then a bounded solution (j) to problem (|^.5p -([X" 
satisfies 

||0!l<T,a + ||V</.||^,a <Ca\\h\\^^a. 

Proof. We already know that 

||0lU,O + ||V^|U,O<C||/l|U,a. 

Then we may write 

'ip{z) = / (j>^{x,z)dx, 
and differentiate twice weakly to get 

ip"{z) — 2 / (j)1{x, z) dx + 2 / (j)zz4>{x, z) dx. 

JR JR 

We have 

/ <j)zAdx= ( (j)ldx+ [ {1 - pwP-'^)(j)^ dx - [ B{(j))(j)+ [ 

JR JR JR JR JR 

Integrating by parts once in x we find 

B((l))(l) = / [-(&i(/))^0^ - {b2(t))^(t)z + b3(t)x(l} + bicjjzcj) + b^cj)^ 

+ 0^) dx. 

Because of the orthogonality conditions (|2.6p we also have that for a certain 7 > 
/ (l)ldx+ I {1 - pwP~^)(l)^ dx > I {(j)l + (j)^)dx. 

JR JR JR 



h(j). (2.10) 



(2.11) 
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Hence, reducing Sq if necessary, we find that for a certain constant C > 

'/'"(•z) > -C I h^ix, z) dx, 

so that 

-^"(z) + J^(z)<-e-2<^l^l|l/.||L. 
4 a 



Since we also know that ip is bounded by: 

mz)\ < -\\h\[i^, 



a 



we can use a barrier of the form '^'^{z) = il/j|/i||^jje + pe^°^, with M sufficiently 

large and p > arbitrary, to get the bound < V ^ ''/'^ for z > and any a < ^ = oq. 
A similar argument can be used for z < 0. Letting p we get then 



/ 



Elliptic estimates yield that for i?o fixed and large 

\<i>(x,z)\ < ae-''l"l||;i||,,, for \x\ < Rq. 

The corresponding estimate in the complementary region can be found by barriers. For 
instance in the quadrant {x > i?o, ^ > 0} we may consider a barrier of the form 



(x,z)^M\\h\Lae-^--+-'^ 



pe'' 



with p > arbitrarily small. Fixing M depending on Rq we find the desired estimate 
for II 011 a letting p ^ 0. Arguing similarly in the remaining quadrants is similar. The 
corresponding bound for ||V(/>||cr,a is then deduced from local elliptic estimates. This 
concludes the proof. |— | 

Notice that for a general right hand side h equation of the form L{(f)) + B{4)) = h 
with the orthogonality conditions imposed as above does not have a solution. On the 
other hand the problem 

L{(t)) + B{(t>) = h + c{z)w.^+d{z)Z, in (2.12) 

under orthogonality conditions 

•{x,z)wx{x)d^.{x) ~ {) ^ I 4){x, z)Z{x)dfi{x), for all z € M. (2-13) 



has a solution in the sense that for given h one can find {(j), c, d) satisfying (|2.12[) - (|2.13p . 

Corollary 2.3 There exist C > 0, r/o > 0, 6a > as in Lemma \2. SX and ao > for 

which, if \\h\\a.a < +oo, a G (0, 1), a G [0, ao), and dfi{x) = p{x) dx is such that 

e""'"'[(k.| + I^DIp.xI + + \Z,\)\P.\] dx < 770, (2.14) 
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then a bounded solution (j) to problem 12^ - ^2. satisfies 

||0IU,a + ||V0|U,a<C||/l|U,a. (2.15) 

Moreover we have 

|c(z)| + |d(z)| <C||/i|U,,e-°l^l. (2.16) 

Proof. To find a priori estimate ()2.15p we have to find bounds for the coefficients c(z) 
and d{z). Testing equation (j2.12p against Wx and integrating with respect to dfi(x) we 
get 



4>zzWa: dfj.{x) + / Lo{(t))uJx dfi{x) + / B {4>)wx dfi{x) = I hwxdfi{x) 

+ c{z) / wld^i{x). 



Let us assume that ||/i||cr,o < +oo and that is a bounded solution. Integrating by 
parts and using Lq{wx) = and the orthogonahty condition (|2.13p we get 

c(z) / wlpdx= / B{(f))wxpdx + / <j>{2wxxPx + WxPxx) dx - / hwxdx. (2.17) 

JR JR JR Jr 

To estimate term J^B((j))wxP,dx we use an argument similar to that of (j2.11[) . and to 
estimate Jj^ (j){wxp)xx dx we use (|2.14p to get 



\c{z)\ < C / \hWx\ + C{6o + V0)i\\yna,0 + ma,o)- 
JR 

From the a priori estimates (Corollarv l2.ip applied to (|2.12p we know that 

||V0||^,O + UWafi < Ci\\h\\a,o + \\c\\oo + \\d\\oo), 

since 

lk.e'^l^^l|U,||^e"l"l|joo <oo, 
for cr e [0, 1). Thus, reducing Sq, rjo if necessary, we find 

||c|ioo < C(|l/i|U,o + ((5o + 77o)NI|oo). 
Testing the equation against Z and using exactly the same argument we find 

ll^lloo <C(||/ii|„,o + ((5o+77o)i|c||oo). 

Hence 

\\dZ\\,,a + WcwxWa.o < C\\h\Uo, 

and the estimate 

||V(/)||,,o + ||(/)||.,o < C||/i||,,o. 

follows. 
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Finally, if we additionally have ||/i||cr.a < +00, we obtain that 

\hW:c\<C\\h\\^^ae-''^'^. 



The same procedure above and the a priori estimates foimd in CoroUarv 12.21 then yield 

\ciz)\ + \d{z)\<C\\h\\,,,e--\'\, 

from where the relation (|2.15|) immediately follows. |^ 

Concerning the existence of bounded solutions of (|2.12p - (|2.13p we have: 

Proposition 2.1 There exists numbers C > 0, Sq > 0, rjo > such that whenever 
bounds l\2.14l hold, then given h with \\h\\„^a < +00, a G (0, 1), a € [0,ao), there 

exists a unique bounded solution (f> — T{h) to problem 12^ - ^2. 13^ which defines a 
bounded linear operator of h in the sense that 

||V0||a,a+ ||0||a,a < C\\h\\„.a. 

Proof. We will first consider solvability of the following problem 

{L + B){(j)) = h, inR2, (2.18) 
in the space V, where ip G V it \\^\\a,o < oo, a G (0, 1] and 

ip{x, z)vuxix)p{x) dx — — / ^p{x, z)Z{x)p{x) dx, for all z e R, (2-19) 

where the density p{x) satisfies the hypothesis of CoroUarv 12.31 We claim that given 
h G V there exists a unique solution (f> of (I2.18P in V. We will argue by approximations. 
Let us replace h by the function h(x, z)x(-r,b) i^) extended 2i?-periodically to the whole 
plane. With this right hand side we can give to the problem (|2.18[) a weak formulation 
in the subspace of H]^ C of functions that are 2i?-periodic in z. To be more 

precise let 

/oo j*R poc 
/ Wtp ■ dzdx + / / ifjrjdzdx. 
'OO J-R J -00 J-R 

By W we will denote the subspace of functions in H}^ that satisfy (|2.19p . Then (|2.18[) 
can be written in the form 

/oo /^R 
/ h^, i^eW, (2.20) 
'OO J-R 

where A : W W is defined by 

/oo pR poo pR 

/ (\/<j>Wij+bi(l)^ij^+b2(l)xi^z)dzdx+ / (1-65)01^ dzrfx, V' G VF, 
00 J-R J -00 J-R 

and K : W W is a linear operator defined by 

/ 00 j*R. poo j*R 

/ [{bix+h2.~h:i)4>.,-hi(t>A^dzdx-p / <t>i^wP-Uzdx,^pGW. 
-00 J-R J -00 J-R 
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Using (|2.4p . (|2.14p . and the fact that Hw^^^^-* || (p_i),o < oo one can show that the 
operator A is invertible and the operator K is compact. 

From Fredhohn alternative and Lemma 12.11 which in addition can be extended pe- 
riodicaUy to a miique solution (j) Cz V, of (j2.18|) with h replaced by h(x, z)x{-R.R)i^)- 
Letting R — > +00 and using the uniform a priori estimates valid for the approximations 
completes the proof of the claim. 

The existence of a solution to (|2.12[) - (|2.13p as well as the rest of the Proposition 
follows from this claim. Indeed, given h such that ||/i||cr,o < 00 by Ilv{h) we will denote 
the orthogonal projection of h onto V (in the sense of LF'{pdx) as indicated by (|2.19p '). 
Using the claim we can solve then the following problem 

{L + Bm=\lv{h). 

Now we only need to chose functions c(z). d{z) such that 

(/ - Tlv)[{L + B){dp)] = {I- Ilv){h) + c{z)w^ + d{z)Z. 

This ends the proof. 
□ 

Remark 2.1 Incidentally, Lemma [2TT] helps us to sketch a proof of existence of Dancer 
solutions (essentially that in [7]) as follows. Consider the bifurcation problem 

AAw-u + w^'^O, inM^ 

with A > 0, which we write, for wxix) = w^^/Xx), as 

AA0 + {pwl'^ - 1)0 + N)^{4>) = 0, 

with N{(f)) = [wx + (j)Y ~ ^P''^\^^4'- We consider the space X of all functions (/) with 
||0||o-,o < +00 which are T = 27r/v'A]"-periodic and even in the z-variable. The operator 
AA — 1 has a bounded inverse in X and thus the problem gets rewritten as 

0+(AA-l)-i(p«;rV + A^A(0)) = 0. 

The derivative of this operator in (/) at A = 1 and (p = is just 

/+(A-l)-i(pu;P-i), 

which has the form: I + K, where K is a, compact operator. Crandall-Rabinowitz 
theorem thus yields that bifurcation takes place at A = 1 since Lemma 12.11 implies that 
this linearization has a simple eigenvalue at A = 1 with eigenfunction Z{x) cos{^/Xiz), 
and a branch of nontrivial solutions in X constituted by a smooth curve S 1-^ {^(S), (j^i^)) 
with (/)(0) = and A(0) = 1, and ds(j>{0) = Z{x) cos{y/Xiz). Scaling out X{S) we obtain 
then solutions wg to the problem for A = 1 with period 2Tr/^/XsX and even in the z- 
variable. With the aid of barriers, we find that these solutions have an expansion of the 
form 

ws{x,z) = w{x) + SZ{x) cos(Va7^) + 0((5^)e~l^l, 
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for all small S. By the standard theory we have that derivatives of this solutions satisfy 
the corresponding relations 



ws,x{x,z) = Wx{x) + SZx{x) cos{^/>^z) + 0{S'^e~^''^), 
ws,xx{x, z) = Wxx{x) + SZxx{x) cos(VAiz) + ©(^^e"!"^!), 

ws,zix,z) = -v/A7(5Z(2;)sin(v/A7z) + 0((52e-l^l), (2.21) 
ws,zz{x,z) = -Ai(5Z(x)cos(Va7z) + 0(5^e~l^l), 
W5,xzix^z) = -\J%Zx{x)smiy^z) + 0{S^e-\''\). 



3 Linear theory for multiple bump-lines 

In this section we will choose a collection of approximations W which will constitute the 
basis of our construction of a solution to Problem p.ip . Let a be a small positive number 
and A: be a fixed positive integer. Let us consider functions fj e C^(E); j = 1, . . . ,k, 
ordered in the sense that 

h{z)<h{z)<---<fk{z), VzGK. (3.1) 

These functions represent approximate locations of the bump lines. Asymptotically 
these lines arc straight lines with slopes proportional to a and whose mutual distances 
are at least 0{dr), where 

= log-. 

a 

Thus, in addition to p.ip . we assume that for certain large, fixed constant A/ > these 
functions satisfy 

fj+i{z)~fj{z)>2d,-M, j = l,...,fc-l, (3.2) 

Here and in what follows we use the following weighted norm for those functions that 
depend on z only 

\\g\\e„a :=||e^°"Wg|U, 

where g : R M™ and 6*0 > is a fixed number to be determined later. Notice that 
the factor 6*0 in front of the exponent a\z\ in the definition of the norm will be taken 
small but fixed independent on a. For convenience we will denote 

f = /2, ■ • ■ , fk)- 

We assume that 

\\i"\\9„a<Ma\ (3.3) 

We will further suppose that the locations of the bump lines are symmetric with respect 
to X axis, which in terms of f means, 

f(z) = f(-z), VzeR. (3.4) 

Let us observe that under these assumptions the numbers 

(3, :=/j(+cx3). 
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are well defined and 

\^j\<Ma. (3.5) 
In fact, from p.3p - p.4p we see that there exist numbers Bj such that wc have 



/,(z) = /3,|^|+i?, + 0(e-««"l^l), 



as z +00. 



From it follows also that: 

pj+i ~ > 0, Bj+i - Bj > 2d, -M, j = 1, . . . , fc, (3.6) 

and 

\\\f-\ - f^Alooc. < Ca, j = l,...,fc. (3.7) 

In addition we will assume that with for some 61, 61 > 9q > with 60 defined as above 
we have 

(3j+i - > e^a. (3.8) 

With functions fj we will associate a change of variables defined as follows: for any 
even function f(z) such that 

\\n\eoa<Ma^ I3^f'{^), 

we set 

X(x,z)^ ^_Ik L=, Z(x, z) ^ \z\^l + {Pvia\z\)y + ' ^^'^^ ^ , 

^l+{(3rj{a\z\)y ^l + {p^{a\z\))' 

(3.9) 

where 77 is a smooth cut-off function such that ri{t) = if t < Ti, ri{t) ~ 1, t > T2, 
where Ti, T2 will be chosen later independent on a. With / = fj and f3 = Pj, satisfying 
(ISH)-(IS3), we will denote X = X^, Z = Z^. 

Let us consider Dancer solutions ws{x, z), even in z, which can be expanded at main 
order as 

ws{x,z) =w{x)+SZ{x)cos{^/%z) + 0{5^)e~'^''^, \d\ < 5a. 
Given /3 G M and S, \S\ < 60 wc define the function 

W5,i3ix, z) = ws{X, Z), whith (X, Z) defined in (3.10) 

With this definition ws^p is an exact solution of the problem in the region where 
r/(a|z|) = 0. If we set /(z) — /3|z| + B then the function 105^^(0;, z) is also an ex- 
act solution of the problem where 77(q;|z|) = 1. Indeed, restricted to the half-plane 
K'^ = {z > 0}, the latter situation corresponds simply to rotating and translating the 
axis of the solution wg from the z axis to the line x = Pz + B. 

We will also consider functions Cj € C^(R), j = 1, . . . ,k, such that 

l|e||eoc + l|e'||0o« + ||e"||eoa < Ma^, (3.11) 

where 

e = (ei, 62, . . . , gfc). 
unctions ej are eve: 

e(z)=e(-z), VzeM. (3.12) 



As in (|3.4p we shall assume that functions ej are even 
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We are now ready to set up our first approximation. Given functions fj, ej satisfying 
dSUl-lim, ((3lT|) - ((37T2)) . and small numbers (Sj, j = 1, . . . , A: we let 

k k 

^{x,z) = ^ws^,fi^{x,z) + ^ej{z) Zj{x,z) , (3.13) 
i=i i=i 

where Pj — fj{+(X)), j = 1, . . . , k, and 

Zjix,z) = Z{Xj). 

In the sequel we assume 

\\S\\<Ma, (5i,...,4). (3.14) 

We want to develop a theory similar to that in the previous section now for the operator 

C{(f>) = A(j> + (p\]P-^ - 1)0. 
More precisely we define the weighted norm 

.. . fe .1 



a\x-fj{z)\-0oa\z\ 



and search for a bounded left inverse for a projected problem for the operator C in the 
space of functions whose || • Ho-^eoa.* norm is finite. Thus we consider the problem 

k 

Ciip) = h + Cj{z)r]j'Wj^x + dj{z)'r]jZj, in R^, (3.15) 

where now we do not assume necessarily orthogonality conditions on 0. Here and in 
what follows we denote 

Wj{x,z) = w{x - fj), (3.16) 

Wj^x{x,z) = w'{x - fj), (3-17) 

Z,{x,z) = Z{x~fj). (3.18) 

We will now define several cut-off functions that will be important in the sequel. To 
this end let ry^(s) be a smooth function with ri^{s) = 1 for |s| < a and = for \s\ > b, 
where < a < 6 < 1. Then, with = log we set 

,b 
,b 
,b 



Pj{x,z) 




Vj{x,z) 




Vj{x,z) 




r]^{x,z) 





25- 


1 


25 




26- 


1 


26 




2'- 


1 


27 




28- 


1 


28 



26- 


1 


26 




27- 


1 


27 




28- 


1 


28 




29- 


1 


29 



(3.19) 



We will prove the following result: 
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Proposition 3.1 There exist positive constants ao, Sq > such that if l\2.4l is satisfied, 
and if a € [0,q;o), <7 G (0,1), and constraints l\3.1]i - ^S.6\i . hold, then problem 

l\3.15^ has a solution 4> = T{h) which defines a linear operator of h with \\h\\cr^eoa,* < 
+00 and satisfies the estimate 

\\4>\\cr.Boa.* + II V0||CT,eoa,* — Cff 1 1 ^1 1 (T,eoa,* • 

In addition function (p satisfies the following orthogonality conditions 

{x, z)wj_x{x, z)pj{x, z) dx = = / (j){x, z)Zj{x, z)pj{x, z) dx, forallz^M.. 

Jr 

(3.20) 

Besides, the coefficients Cj{z) and dj{z) in H3.15\i can be estimated as 

k 

^(|c,(z)| + |d,(z)|) < C||;i||.,,„„,.e-«»"l^l. (3.21) 
j=i 

Proof. The main idea in the proof of this proposition is to decompose problem (|3.15p 
into interior problems that can be handled with the help of the theory developed in the 
previous section and an ea;<enor problem and then glue the solutions of the subproblems. 
Form the definition of the functions r]j , rf^ we have 

^3^J = ^7' ^Ph = ^j- (3-22) 
We search for a solution of (|3.15p of the form 

k 

Substituting this expression into equation (|3.15[) and arranging terms we find 
fc 

^77,[A0, + (1 -pWP-i)0, - Cj{z)wj,x - dj{z)Zj - h] 

k k 

+ [a^j + (1 - pWP-i)V> -{^-Y. ^0 ^ ~ + Ar;j0j) 

We will denote 

hj=vpi, =p77+(wf i-WP-i). 
Let us observe that in the support of rjj we have, using p.22p . 

Then we find a solution to problem p.lSp if we solve the following linear system of 
equations 

— (1 - pwp^)(j>j + rj(j)j = hj — p\}P~^T]~-tJj + Cj{z)wj^x + dj{z)Zj, (3.23) 
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in K^, for j = 1, . . . , fc, and 

k 



1 - E'?^)"'"']^ = (l - E'?^'^ - E(2Vr?,A70, + Arj.cl,,), (3.24) 



in M^. To solve equations p.23p we denote <j>j 
write the equation for (f)j 



'qlii and use ^^-^^ to 



- (1 - pw^-y^, + rj0, = + - WJ'-i) + r, + p(l - E 

+ CjWj,x + djZj. 

(3.25) 

We observe that equation (|3.24p written in terms of 4>j has form 

- [1 - E'?^) = (1 - E^j)^ - E(2V??j^0, + A77,^,), (3.26) 

since for instance \/rij'\7{j]Jrlj) = 0. 
Let us denote 

L,(,^) = A0-(l-p«;f 1)0, 
and consider first the auxiUary problem 



Lj{(j)) ^ h + c{z)wj^x + d{z)Zj, in R^, 
under orthogonality conditions 

<{x, z)wj^x{x, z)pj{x, z) dx ^ ^ / <j>{x, z)Zj{x, z)pj{x, z) dx, for all z e 

Jr 

For future references observe that if Pj is sufficiently small then 



(3.27) 



(3.28) 



(3.29) 



We want to solve p.27p ~ (|3.28p using Proposition 12. II To this end we consider the 
natural change of coordinates 



= X-fj, Zh->Z. 



and set 



(j){x,z) = 0(X,Z). 

Direct computation then shows that problem (|3.27p - (|3.28p is equivalent to 
L{4>) + Bj{4>) = h + c(z)wx(x) + d{z)Z{x), in 



(3.30) 



under orthogonality conditions 

/ 4>{x, z)wx(x)p(x) dx^Q 



(x, z)Z(x)p(x) dx, for aU z e M, (3.31) 
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24 



where 



X \ 2^ - 1 , 2^-1 



26 
Here 

and 



- (I)'*" l-^^W' l5-^"w-<"^> 



The operator Bj satisfies the assumptions of Proposition l2.f 1 since from p.3p we have 
using the notation of (|2.4[1 and denoting the vector of the coefficients of Bj by hj : 

||b,|| <C(||/j||o,, + ||/;||oo) <Ca. (3.33) 

Problem (|3.30p - (|3.3f p has a unique bounded solution (j> = Tj{(j>) where Tj is the linear 
operator T predicted by the proposition for B = Bj. Besides, if ||/i||o-,q < +oo then we 
have the estimate 

going back to the original variables we see then that there is a unique solution to (j3.27p - 
(|3.28p . 4> = Tj{h) where Tj is a linear operator. In addition we have the estimate 

where 

Moreover, the coefficients c and d are estimated using (|2.I6p by 

|c(z)| + \d{z)\ < C||/i||,,(,oo,je-««"l^l. (3.34) 

Let us observe that given ijj we can recast the equations p.23p for (j)j as a system of 
the form 

j ^l,...,k, 

(3.35) 

We will solve next equation (|3.26p for -0 as a linear operator 

where $ denotes the /c-tuple <& = {(f>i, . ..,4>k)- To this end let us consider first the 
problem 

A^P - {I ~ 9)i' = g inK^. (3.36) 

where 

k 
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Observe that if the number is large enough then 9 is uniformly small, indeed 9 = o(l) 
as a — > 0. Let us assume that g satisfies 



\gix,z)\ < yi^e-'^l^-/^(^)l-««"l^l, 

for some < fi < 1. Then, given that if rf* is sufficiently large then number 9 is small, 
and also that (|3.3|) holds, the use of barriers and elliptic estimates proves that this 
problem has a unique bounded solution with 



Thus if we take 



we clearly have that 



\g{x,z)\ <C ||/l||a,eoa,* {\\4>j\\a,eoa,j + \\^<i>j\\a.eoa.j) 

k 
J = l 

and hence equation p.26p has a unique bounded solution 

which defines a linear operator in its argument and satisfies the estimate 



oc* + O{l)^{\\(j)j\\a,eoa,j + l|V^j|U,eoaj) 

i=i 



X 



-;'|2:-/3(z)|-9oq|z| 



In addition, we find that 



|^||a,eoa,* +o(l)^ (II 

+ l|V0j|U.0„aj) 

i=i 



(3.37) 



(3.38) 



Now we have the ingredients to solve the full system (|3.25p - (|3.26p . Accordingly to 
(|3.35p we obtain a solution if we solve the system in $ 



+ Tj (^7*(<i>, 0)x, + r,4>j) = T, {h, - ^-^{0, h)xj) , J = 1, . . . , fc, (3.39) 
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where 



Xj 



We consider this system defined in the space X of all functions <& such that the 
norm 

k 

is finite. System p.39p can be written as 

where A and B arc linear operators. Thanks to the estimates for the operators Tj and 
the bound p.37p we see that 

\\B{h)\\x < C\\h\Ueoc..*- 
On the other hand we have that 

k 

Umix < c[Y,\\vJ^{<^,0)xj\\a,eoa,j+Y.\\r,^j\Ueoc..j ■ 
Using estimates p.37p and (|3.38p we find 

k 

Y.hJ'^{^,0)XjLM < oiimWx- (3.40) 

J=l 

From the definition of Vj and (|3.14p we get Ht-j Hoo = o(l) which implies 

k 

Y.\\^j^3\\-M<o{i)\m\x. 

Summarizing the last estimates we obtain 

MWILy <o(i)||a>||x, 

hence the operator ^ is a uniformly small operator in the norm || • \\x provided that a 
is sufficiently small. We conclude that system p.39p has a unique solution $ = ^{h), 
which in addition is a linear operator of h such that 

mh)\\x <C\\h\\,,eoc.,*- 
Thus wc get a solution to problem p.lSp by setting 

k 

= E + (l - E h). (3.41) 

Using p.29p we get pjcfi = pj(f>j hence from and (|3.28p we obtain (|3.20p . Estimate (|3.2ip 
follows directly from (|3.34p . The proof of the proposition is complete, q 
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4 The nonlinear projected problem 

Let us recall that our goal is to find a solution of the problem 

S[u] -.^ Au + uP - u ^ Q, p>2, inM^ 

which is close to the function W defined in (j3.13p . We will denote by E the error of 
approximation at W: 

E := 5[W] = Aw + W** - W. 

A main observation we make is that, under the assumptions (|3.ip - p.6[) . (|3.14p . W is such 
that 

^*:=ll^[W]||a,eoa,*<C«'~'", (4.1) 

where 6*0 is defined in (|3.8p . We will postpone the proof of (|4.ip for now and for the 
purpose of the present section we will simply accept it as a fact. We look for a solution 

to our problem in the form 

where (/> is a small perturbation of W. Thus the equation for u is equivalent to 

£(0) = £:-iV(0), inR2, (4.2) 

where 

£((/.) :=A0+(pWP-i- 1)0, (4.3) 

and 

iV(0) = (W + 0)P - WP - p V- (4.4) 



Rather than solving problem (|4.2p directly we consider an intermediate projected version 
of it, 

fe 

£((/)) = E- N((l)) + Cj{z)iywj^a: + dj{z)T]jZj, in (4.5) 

where rjj are the cut-off functions defined by (|3.19p for problem (|3.15p . 

We will establish next that this nonlinear problem is solvable with similar estimates 
to those obtained in Proposition in l3.1l for problem (|3.15p with h — E, namely we show 
the following result. 

Proposition 4.1 There exist positive numbers ao,So, ctq, such that for any number 
<J G (0, (To), a G [0,q;o) and any f, e and S satisfying constraints jp - (|g.6'p . I\3.14l , 
Problem l \4-5^ has a solution (j) with \\(l)\\a-,eoa,* < Ca-'a^~^'^ , where a' G (ct, 3(t/2) such 
that 



{x, z)wj^xix, z)pj{x, z) dx ~ ~ / (j){x, z)Zj{x, z)pj{x, z) dx, for all z 

Jr 

(4.6) 

The coefficients Cj and dj can be estimated as follows: 

k 

J2i\^,{z)\ + \d,{z)\) < C.,«2-2-'e-^°"l^l. (4.7) 
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Proof. Let us observe that we obtain a solution of the problem (|4.5p if we solve the 
fixed point problem 

= r(i?-iV(0)) (4.8) 

where T is the operator found in Proposition l3.1l Assume that H^jllo-.eoa,* < I7 J = 1) 2- 
We have that 

I N{cl>,) - 7V(0i) I < + |(/.2|) 102 - 0l|, 

and hence 

II iV(02) - N{cj)i) |U,e„a,* < C( UlWaA^c* + Whh.Boa.*) ||02 - (/-l |U,e„a,* , (4.9) 

in particular 

II A^(<^) iu,eoa,* < mio,,^,,- 

Then, from (|4.1[) . the following holds: for each fixed a' G {a, 3a/2) there exists a number 
V > such that for all small a the operator is a contraction mapping in a region of 
the form 

^ = {0/ \ma.e„c.,*<i^a''^"'}, 

and hence a solution of the fixed point problem (|4.8p in B exists. Furthermore (f> solves 
(|4.5[) . and by (|3.20p we find that satisfies (|4.6p . The proof of the proposition is 
complete, q 



5 Estimates of the error of the initial approximation 

In this section wc will show estimate (j4.ip announced above for the error 

E = 5[W] = Aw- W + W^. 

We will denote 

Ei = AV, E2 = -V + \}P, (5.1) 

and let for j = 1, . . . , fc, 

,Jl + {l3Ma\z\)y ^l + {/3M^\z\)y 

With a j fixed we will set for brevity of the notation fj = f, Sj = 6, (3j = (3, Xj = X. 
Zj = Z and let 

ws,pix - f{z), z) ws{X{x, z), Z{x, z)). (5.2) 

We observe that X{x,z) ~ X(.t, — z) and Z{x,z) = Z{x,—z) and therefore we only 
need to consider the error of the approximation assuming z > 0. For brevity we will 
denote 

1 1 



ao = >/l + iPv) 



2 



ai = — = 



a2 — o-i — n , cia — ^2 
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Using this notation wc can write, whenever z > 

X = (cc — /)ai, Z — zao + PtjX. 

We have that 

= ai, X^ = -/'fli + Xa2, 

Zj. = Prjai, Zz ~ ao + api^zrj'rjai + a/Jyy'X + /JyyXz 

= Oq + a(3'^ zrj'rjai — /'/?7/ai + ajSrj'j. + jSrja^T^^ , 

and therefore 

An := X2 + x2 = a2(l + (/')") - 2/'aia2X + 0^x2 

= 1 + al{{f'f - /3V) + 0(a3)(l + X2)e-^''"l^l, 
^12 = A21 X^Z^ + X^Zz = Pi^al - fai{ao + a/Sij'X + IStjX,) + aaXZ^ 

A22 -.= 11 + 21 = 1 + 2/?77(/377 - /') + ((/')' - + 2aP^zT]'r] + 2aoaPr]'x 

+ 0{a^){l + X^)e~^°"^'^. 

Similarly, wc get 

h := Ax = -f'ai + a2(-/' + X^) + Xag, 

I2 ■■= AZ = 2(3'^mj7j'ai + 13"^ az{ri' r^ai)' + PijAX + 2a[3Xz'n' + 

Denoting A = {Aij — 6j), I = {h, I2) we have for any function u = u(X(a;, z), Z{x, y)), 

l^x,zU = Ax,zu + Tr[A ■ Hessx,z(w)) + 1 ■ Vx,zu, 

From (|3.3p - p.7p wc also have the relations 



A = 

and 



0(a2) + 0(a3)(l + x2) Oia) + 0{a^)il + x") 
0(a) + 0(a3)(l + x4) 0(a2)(l + x4) 



1 = {0{a^) + 0{a^)X,0{a^) + 0(a3)x)e -Soal^l , 
Using now the fact that 

k k 

W(x, z) = Y^ (Xj, Z,) + ^ e, (z)Z(X,), 
for each j = 1, . . . , fc and (|2.2ip and (|3.1ip we obtain 

k 

Aw = ^ Ax^.z.u'i, (Xj, Zj) + E12, (5.3) 

where 
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We now turn to computing £'2- For brevity we will set: f{u) ~ —u + , p > 2. We fix 
a j, and li \ < j < k we define sets 



= (x, z) 



fj-i{ez) + fj{ez) fj+i{ez) + fj{ez) 



while when j = 1 or j = k we set 



^1 = 1 ix,z) 
Ak = \ ix,z) 



-00 < X < 



/2M + /1M 



fk-l{£z) + fk{£z) 



< X < 00 > . 



For X € Aj, with I < j < k, fixed we write 



/' (u'a, ) (W - ,;3, ) - XI / ) 



+ 2/"K-,ft)(w-^i'5,,ft.)' + £;2j. 

We first assume that p > 2. Since p > 2 we have 



\E2j\<c[ 



maxe 2 



< Ca(P"^^e"5(p-3)eoa|z| j^^g-<T|x-/.|g-i(3-<T)|/,-/.| ^ ^3 ^ g-3eoa|^|g-3b-/.| 



k 



1=1 



Thus when < p — 2 then 



\\E2j\\aA>o.,* < CaP-\ 
Going back to estimation of E2, we have 

/' ( ) (W - ^5, ) - E / (^-^-^^ ) 



= maxO(e"(P-i'l^-^^l-l^--^'l) +maxO(e-Pl^-^-l) 



+ O(a2)e-eo«klg-|x-/,l^g- 



-(H-/i)|a;-/i| 



i=l 
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When < a < p — 2 we have for i j 

(p -l)\x- J,\ + \x- J,\ = a\x -,f,\ + {p-l- a)\x -f,\ + \x- h\ 

> a\x - .f,\ + min{l, (p - 1 - a)}(|.T -f,\ + \x- /,|) 
>a\x^f,\ + \f,-M 

> a\x - /,| + (lA - Pj\\z\ + 2d, - M), 



and 



p\x - .h\ > a\x - .M + - + 2d, - M). 



Using IPi ~ l3j \ > 9oa wc get for x £ Aj and p > 2 

AW + /(W) = 0(a2-2'^)e-''l^-^'(^)l-^°"l^l +£^2^, (5.4) 
hence in Aj we have 

\\E\U0„a,* < Ca''-^r (5.5) 

When p ~ 2, we triviaUy have 

= \\E2j\\aA,c.,* <CaP-\ 

and 

\x - /jl + N - h\ - ct\x - f,\ + (1 - a)|a. - /,| + |x - h\ 
><j\x-fj\ + {l-a)\J,-n\ 
> cr\x - /j| + (1 - (T){eia\z\ + 2d, - M). 

Since 9i > 9q, we may choose a smah such that (1 — (7)61 > Oq and hence (|5.5p also 
holds. The ends the proof of the estimate, q 

6 Dependence of the solution on the parameters 

Now we will study the dependence of the function </) found in the proposition above 
on the parameters f , e and S. More specifically we are interested in establishing the 
Lipschitz character of as a function of variables (f",e"). We will begin with the 
observation that the error term 5[W] can be written as follows: 

5[W] = Ei{i", e"; f ; 6) + E2{i' , e', f , e; S), 

where for functions g, h which are even and a vector parameter 5 ~ {Si, . . . , Sk) such 
that 

||g"|U,||h"||„ <Ca^ ||^||<Ca, (6.1) 

we have defined: 



Ei{g", h"; f; 5) = - ^ g^a^l{l,f3,) ■ Vx.,z.?«5.,ft + ^ h"Z{X,), 

^=l ^ ^=l ^ ^g2) 

E2{i', e'; f , e; 6) = 5[W] + ^ n'a,i{l, A) • Vx.,z.u.5.,^^ - ^ e"Z(X,). 
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Notice that E2 depends only on f , f, e, e', S and it has exactly the same expression as in 
the previous section. Given g, h, f , e, f , e, 6 let us now consider the following problem 

k 

£(</)) = £;i(g", h"; f; d) + E2{i', e'; f , e; S) - iV((/)) + ^ Cj{z)T^jWj,, + dj{z)TjjZj, in . 

(6.3) 

The solution of this problem can be obtained by the argument of Section 21 In fact the 
analog of Proposition WA\ yields H^Hq.o-,* < Ca^^*^ , where cr' e (cr, 3<t/2). 

We will consider functions gi, h.^, f^, e^, f^, e^, and vectors Si, i = 1,2 such that 

\\Si-S2\\eoa+a\\i[-t2\\e„a+a^\\ii-h\\^ < Ca\ 
||h'/ - + ||e'i - e'alleoa + ||ei - e2||e„a < Co?, (6.4) 

II (5i - (5*211 < Ca. 

First we want to show that functions Ei and E2 are Lipschitz functions of its variables. 
To make it more precise we will distinguish the norms taken with respect to different 
functions f . Thus we will denote in this section: 



||</'IU..ca,f :=|(5]e-'^l--/^(^)l-«o.|^l] cj, ^, f = (A,. ..,/,). 
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Lemma 6.1 Under the assumptions {6.4^ we have the following estimates: 

||i;i(g'/, hi; fi; 5i) - i;i(g^', hs; f2; l2)IU,eoa,fi 

< C(||g'/ - g^'IU + ||h'/ - h'%„^) + C(max ||gf ||,„„ + max \\h':\\e,^)\\t^ - t2\\oo 

i—1,2 'i— 1,2 

+ Cmax||gnieoa||^i-52||, 

1—1,2 

(6.5) 

||£'2(f{, e'l; fi, ei; (5i) - -B2(f2, §2; f2, 62; (52)|U,eoa,fi 

< C(ai-2-||f{ - f^||e„„ + lle'i - ey|e„„) + C{a'-^''\\h - h\\oo + ||ei - e2|U„a) 

+ Ca^-^''\\A-52\\- 

(6.6) 

Proof. First notice that to prove Lipschitz estimates (|6.5p - (|6.6p it suffices to show 
corresponding inequalities varying one component at a time (for instance considering 
gi, i = 1,2 and fixing the rest of the parameters). When f and f are fixed then the 
estimates are rather standard. We will therefore concentrate on the case when all 
parameters except f and f are fixed. 

We will indicate the dependence of the change of variables Xj , Zj on f by writing 
Xj(f),Zj(f). Similarly, we will write Aj{{,i') for the components of the matrix Aj. 
Finally we will set 

l,(f,f')=lj(f,f',g")+5>«i(l,A). 
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We will gather first some estimates that are important in the proof. We have 

HX.ih) - w'iX,i{2)\ + h(X,(fi) - w'iX,{{2)\ < C\wiX,if^)\\h - f2|, 

|l,(fi,f') -l,(f2,f')l < Ca^lfi -f2|, (6.7) 

|l,(f,f0-l.(f,f2)l <C^«|f(-f2l- 



From the definition of Ei and (|6.7p one gets (|6.5p . Estimate (|6.6p is somewhat 
tedious but rather standard. For instance let us consider a typical term coming from 
AW: 

\Tr{A. Hess^^ ,z, {ws^ )) (fi X) - Tr (A, • Hess^^ ,z, {ws^ ,0,)) (fs M 
< \Tr{A, ■ i/essx„z,(u;,,.,ft.))(fi,fO - Tr{A, ■ Hessy,^,z^iws^,i,^)){hX)\ 

+ \Tr{A, ■ i7essx,,z,(«^5,,ft))(f2,f^) - Tr{A, ■ Hess^^,,^{ws,^p,)){hM 
= B^+B2. 



Using the (|6.7[) we can estimate the first term above by 

5i <C|A,,n(fi,f0-A,-n(f2,f0lh(Xj)(fi)l 

+ C|A,,n(f2, fOIK(X,)(fi) - ^"(X,)(f2)l 
+ C<5,|A,(fi,fO-A,(f2,fOII^(X,)(fi)l 
+ C5,|A,(f2,fOII^(X,)(fi)-^(X,)(f2)| 
<Ca2|fi_f2|e-«o"l-l|i^(x,(fi))|". 

In an analogous way we can estimate B2 ■ Another typical term we have to deal with is 
(considered in the set Aj defined in the previous section) the following 

\f{w,^^p^){^ - Ws^,p,){h) - nws,^p,){^ - WS^ ,^ J (f2 ) I 

< I /' im, ) (fl ) - ./' {WS, ) (f2 ) 1 1 (W - WS^ ) (fi ) I 

+ l/'(«^^„ft)(f2)ll(W-li'5,,/3,)(fl) - (W- W,^,^J(f2)| 

< Ca^— e-^«"l^l|fi -f2||i«(Xj(fi))|" + Cmaxe-l-^«--^"l|fi -f2||i(;(x,(fi))|" 

< Ca^—e-^^^l^llfi - f2||«;(Xj(fi))|". 

Other terms in the definition of E2 are handled similarly and from this we get (|6.6p . 
We leave the details to the reader, q 

In what follows we will emphasize the dependence of 4> and Cj , dj on parameters by 
writing 

</>W =0(gf,hf;fl,e^;f„e,;j:), 



„(«) _ Ji) f„ll U". Cl X' 

i ~ i 1 '■'■I 1 '-1 1 '^n '■n Oi ) 



cr =cy^(gr,h^';fl,e^;^,e,;<5,), (6.8) 
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Proposition 6.1 Let (f), be the solution of i6.3\) . Then for j = 1, . . . ,k functions 4>, Cj, 
dj are continuous with respect to the parameters g", h",f', e', f, e. Moreover assuming 
j6'.4P ; in the notation of I16.8\ } we have the following estimates 

<Ca-2'^'(||g'/-g^'|U„a + ||h'/-h^'||,„«) 

+ Ca-'-'{a\\i[ - f^||,„„ + lle'i - e',\\g„^) (6.9) 
+ Ca-2"'(a2||fi - falloo + ||ei - eajls^a) 



Ca2-2-'||5i -(5*21 



where a' G (cr, 3cr/2). 



Proof. First observe that it is sufficient to prove Lipschitz dependence considering 
dependence of cf) on different components talcen separateiy. For instance iet us fix 
f, e, f and e,d and let vary g^, hj. Tlien tlie proof of the result follows directly 
from Proposition (|3.ip and Proposition (|4.ip applied to the equation for the difference 
(/)(g'/, h"; f', e'; f , e; S) — 0(g2 , ; f', e'; f , e; S). We leave the details to the reader. 

We will now consider Lipschitz dependence of (f> on f',e', e. Let us fix g and h as 
above as well f and S. We will denote by (/j*^*^ the solution of the following problem in 



k 

£(0(')) = E,{g", h"; f ; S) + E^ifi e',;{, e,; S) - 7V(0(O) + ^ cf\z)^,w,.^^ + df\z)iyZ, 



3 ■ 



It follows from Lemma 16.11 that 

\\E2 (f( , e'l ; f , ei ; ^) - ^2 (f^ , ; f , 62 ; S)\\,M^ 

< Ca'^"^ {a\\{[ - fslleoQ + l|ei - e2||0„a) + l|e'i - ejHeoa). 

Applying now the theory developed for the linear problem for the difference 
we get 

< Ca''^" {a\\{[ - i2\\eoa + l|ei - e2||0na) + \\e[ - ejHeoa)- 



(6.10) 
(6.11) 



,(2) 



(6.12) 



Finally, it remains to consider the Lipschitz dependence of solutions on e and f . This 
case is somewhat more complicated since the norms involved depend on the variables 
Xj , which in turn depend on the functions fj . 

Now, by (Z)^*-* we will denote solutions of 



£f. (0W) = E,{g", h"; f,; S) + i?2(f', e'; f„ e; S) - 7V(0«) 
+ Y,(fiz)7jfw^ + d^{z)rjfzf\ in] 



(6.13) 



with fixed functions g, h, f , e, and function , i = 1, 2 such that 

||fi-f2||oo <C. (6.14) 
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The equation for the difference (p^^^ ~ (f)^'^'^ = (j) can be written in the form: 

£f, = i?i (g" , h" ; f 1 ; ^) - i?i (g" , h" ; f2 ; ^) 

+ E2{i', e'; fi, e; S) - S2(f', e'; £3, e; 5) 

+ p[(w(i))P-i - (w(2))P-i]0(2) + [iV(0(i)) - Ar(0(2))] 

Observe that if we denote f = fi — £2 then we have ||r||e(,c< < C and 

k 



|w(i)-w(2)|<qr|^e 



It follows 

llW«-W(2)|l,,^,f^ <C|lr|U. (6.16) 

Similar estimates hold for other terms involved in the right hand side of (|6.15p . Now, 
we will further decompose 

k 

where 

The Lipschitz character of functions 4'j, j ~ 1, . . . , /c can be established directUy from 
the definition. On the other hand function 4>q satisfies an equation similar to (|6.15p 
together with the orthogonality conditions. This, estimates in Proposition 14.11 and 
Lemma |6 . II vield finally 

110^'' - 0^'^IU,eoa,fi + ||V0(i) - V^^''^Uo,o.,f, < Ca2-2-'||f||^e"''°"l^l. (6.17) 

7 Projection of the error on the elements of the ker- 
nel 

Let us now compute the projections of the error on the elements of the approximate 
kernel. Thus for each j = I, . . . , k we will compute 



n/, := [ SMw4x,)dx, 

:= [ SMZiX,)dx, 
Jr 
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For a fixed j we begin with H/^. . We will use decomposition of 5[W] similar to the one 
in section O 

k 

SM = [^^(A' • Hessx„z,iw5,j3.) + h ■ Vx„z 



i=l 

k 



(7.1) 



i=i 



S 



52j[e,Z(Xi),...,Z(Xfe)] + 53,[W]. 



Using (|2.2ip we obtain: 



w"(X,) 




5r 



Z"(X,)cos(VXrzO -x/XrZ'(X,)sin(^ZO 
VArZ'(X,) sin(V^Z,) -AiZ(Xi) cos(VArZi) 

From the formulas for X.; , Z; and their derivatives we also get 

aim' - Piv') - n + «fi - PW) 



A,: = 





2aoaP,T]% 



0(a3) 0(a3) 
0(a3) 0(a3) 



(1+X2 + X4)e-««"I-I 



= Ajo + Ail + Ai2. 
Similarly we have 

Vx.,z.u;5.A = K(X.),0) + ,5,(Z'(X,)cos(yA^z0,-yArz(x0sin(yA7z0)+O(52)e-l^'l, 
and 



1, 



We observe that the entries of A.; and 1^ are polynomials with respect to X^ and they 
depend on only through Xj. In the sequel we will use also expansion 



cos(vAiZj) = cos(vAiajo2) + 0{a)\Xi\, sin(vAiZj) = sin(vAiaioz) + 0(a)|X,;|. 

Since functions w" , Z and Z" are even functions of its arguments and w' is an odd 
function therefore, changing variables x t—f Xj^ we get 

rr(A, • Hess{ws^^p^))w'iX^) dx = -2^,ajoS, - /j + a%{if^)^ ~ p^r^^)] 

X sin(-\/A]"aio2^) / w' {s)Z' [s) ds 



+ <5,0(<5,+/3,)(/j-/3,77) 
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When i =/= j we have 

with similar formulas for the integrals involving Z'{Xj) and Z"(Xj). It follows then 
We have also 



+ 0(«3)(/; + l)e-«o"l-l^ 

and similarly for i 7^ j we get 

1, •Vx.,z.^^5.,ft^'(X,)rfa^ = -/rO(e-l/^-/'l) + 0(«3)(/; + l)e-«°"l^l. 



To compute the corresponding projections of S2j we notice that 

Ax,,(e,(z)Z(X,)) = e'^{z)Z{X,) + 2e^(2)Z'(X,)X,, + e,(z)A,,,Z(xO 
= e'IZ{X,) + [2e^(-/;a,i + X.a.s) - e,/>,i 
+ e,(/; + X0e-'«"l"IO(a^)]Z'(X,) 

+ e.[l + a?i((/;)^ - Phf) + 0(«')(1 + X? + x4)e-''o"l^l]z"(X,). 

(7.2) 



Then we get 



S2jw'{Xj) dx = -(2e' /' + e^/") / Z'(s)i(;'(s) 



5]{er-(2e,/^aa+e,/;V) 

+ e,[l + afi((/;)2-/3,V)]}0(e-l/-/.l) 
^(e^ + ej;)0(a3)e-''''"l^l. 



Using the sets Ai, I = 1, . . . , k introduced above the projection of S'3_,[W] on w'{Xj) can 
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be written as follows 



/ Ss,Mw'{Xj)dx^ [ \f{w5^^p^)y^Ws,,p,-y^f{ws,J3,)]w'{Xj)dx 



1^3 
k 



I /'(^5<.A)I]e,Z(x0u;'(Xj)dx 
1=1 •'^i 1=1 

I 0{i\i-ws,j,f)w'{X,)dx 

;_i J A, 



1=1 ''^^ 



For 5*3 J 1 we get 



A, 



For a fixed i ^ j, using formulas (|2.2ip . we have 

[pw{Xj)P-''^w{Xi) - wP{Xi)]w'{Xj) dx 

i 

+ Sj f [0(e-(i+^)l'^l)u;(X,) + 0(e-(i+^)l''l)w(Xj)]rfa; 
Jaj 

To compute In let us observe that for any A G M and a < 6 we have: 

pvjP ^{s)w'{s)w{s ~ X) — / w''(s — A)w'(s) = pw^(s)ix'(s — A) L 



[wP{s)w'{s - A) + wP{s - X)w'is)] 
= [pwP{s)w{s - A) + w'{s)w'{s - A) - w{s)w{s ~ A)] , 
hence, changing variables, we get, with sj ~ ^^-^-^^ tj = , Ay = /,; — /j, 

[pwP(s)w(s - Ay) + w'{s)w\s - Ay) - i(;(s)w(s - A^- )] 
+ 0(/3f +/3f)e-l/^-/-l. 
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Similarly, 

/2, = 0(<5,)e-l/^-^'l. 

Using this formula with i ~ 1, . . . , fc and noting that for a fixed j the dominating terms 
in Ssji come from i = j — I and i = j + 1 wc get 

33,1= w[ ^ ) +w[ ^ ) -w[ ) -w[ ^ ) 

i¥=3 

where we have denoted 

Cp = \ixn{e'[w{s/2f + w'{s/2Y]}. (7.3) 

s — *oo 

The remaining term 5*3^2 can be estimated in a similar way. Observing that when 
I ^ j and X € ^/ we have w'{Xj) = iy'(X;)0(e"l-^J"-'^'l), we get that 



S3j2 = 0{a) e 



To compute S^jj, we observe that 



f'{s)Z{s)w'{s)ds = 0, 



hence 



0{e.j5.j) + "Ye,0{e-\f^-f'\). 



Finally wc have 



S,,,^0{e]) + Y.e^,0{e-\f^-f'\). 
Summarizing wc get the following expression 



II/j- = ^//[co + ciGj + 0{5j)] + 2V XittjoSjif- - Pj-q) sin(v XiUioz) 
- 2ci/je;. + 7,(/j - Pfl) + 0{e]) + e,0{5,) 

+ ^(e, + 0(a))e-l-^^-^'l, 



(7.4) 
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where functions "fj, puj and quj are smooth functions of their arguments satisfying 
estimates 



7j = 7,(/;, /3j, 5„ a, z) = 0(<5f ) + 0{a^), 

- qHj{fjJ^,e,, a, S, z) - 0{5,) + O(e-l^-^'l) + 0(«3g-''''"l^l), 



and Cp, Co, ci are constants defined by 



C„= hm e*(u;(t/2)2 +w'(i/2)2), cq = / ^'(sj^ds 



Cl 



t(;'(s)Z'(s)ds. 



Now we compute projection of the error on Z(JLj) denoted by W^j above. We will use 
the expression for S'[W], (|7.1[) . Wc will denote the entries of the matrix Aoi by ^oi,im- 
First we observe that 



Tr(A^ • i?essx,,Zi('u;5,,/3j)Z(Xi) ^Oi,iittiO / w"{s)Z{s)ds 



Si Aoi^ii cos{\/ XiQioz) / Z"{s)Z{s)ds 
L Jr 

2Aoi,i2aiQ\/%_sm{^/)^aiaz) / (Z'(s))^ds 

AiAoi.22aio cos(VAiaio2;) / Z^(s)ds 

Jr 

+ S,0{\\n-PM\eoa+a')e-'<"'\^\ 



Oi- 



(7.5) 



It is convenient to write 
Notice that 



||1,|U + a^'lll'JU < + II/; - A^IU + |X,:|). 



Then we have 

li • Vxi,z.W5i,/3,^(Xj) = (5i/"a,a/3i77sm(V^aio2) / Z"^ + 0{a) 

+ 0(a3)(l + ||/,'-/?,r;|U + <5,)e-^«"l^l. 

When i j we get 

/ rr(A, • iJessx„z,K.,/3j)^(Xj) 

= 0(||/,'-A??l|a+«')e-'«"l^l(l + <5,)e-l^--^^l, 

and 

/ l.-Vx.,z.^i'5.,ft^(X,) = 0(||/riU + ||/;-A??|U + a2)e-^°"l^l(l + <5,)e-l^--^^l. 
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Now we compute the projection of S'2j[W]. We notice first that 
Jr 

where 

do= / Z^{s)ds, di^ I sZ'{s)Z{s)ds, ^2 = / Z"{s)Z{s)ds. 
Jr Jr Jr 

On the other hand 



V / A,^,{e,Z{X,))ZiXj)dx 



+ ^ e, [1 + /raa + ((/D^ - /3,V)] 0(e 

+ ^e,a3(o(|/,^_/^.|) + /;)0(e-«oa|.|)o(g(i+M)l/,-/,l). 

To calculate the projection of S^j [W] we observe that 

/ ^3,[W]Z(X,)dx = ^ / /'(u;,„ft)^e,Z(x0Z(X,)d.T + ^O(e-l/-/^l). 
We have 



On the other hand when I ^ j then 

e.J'{ws,,p,)Z{X,)Z{Xj)dx = 



e 



Summarizing we get 



Ilgj = (e" + Xe-i)do + e-jgoj + e'^gij + ^(eiroy + e-riy + e-V2y ) + /ij, (7.6) 
where functions gij , ruj are smooth functions of their arguments that satisfy 
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Function hj satisfies 



hj — hoj + S J ajiPiT] siii{\/Xiajoz) + 0{a) 

+ 0{{a^l + \\n-l3M)i'^ + Sj) 



with hoj defined in (|7.5p . 

8 Derivation of the reduced problem 

In this section we will derive conditions for f , e and 5 which imply that 

c,(z) = 0, d,{z) = 0, j = l,...,k, (8.1) 
and consequently lead to a solution of our original problem. We recall that with 

we have 

{x, z)wj^x{x, z)pj{x, z) dx = ^ / (j){x, z)Zj{x, z)pj{x, z) dx, for all z e M, 

Jr 

where cf) is the solution of 

fe 

=E- N{(P) + ^ cj{z)7j,wj,x + dj(z)ri,Z, in M^, (8.2) 

(see (|43)) . 

To derive the reduced problem we will multiply (|8.2p by PjWj^x and pjZj, j — 1, . . . , fc 
and integrate over R with respect to x. To begin with let us observe that 



wl^rjjPj dx = 51 / (w')' + 0(a), / ZjryjPj dx = / + 0{a), 

i ' JR JR JR 

hence for (|8.2p to be satisfied we need ioi j ~ 1, . . . , k: 

[E - N{cb) - C{cb)]wj^xPj dx = 0, (8.3) 
[E-N{(l))-Ci(j))]ZjPjdx = 0, (8.4) 



where (j) is the solution found in Proposition 14. II These equations can be written also 
in the form 



n/, = / [N{(f>) + C{(f>)]wj,xPj dx + / E[w'iXj) - Wj^xP,] dx = Q,{f, e), (8.5) 

Jts. Jr 

= / [7V(0) + CmZ.pj dx+ f E[Z{X,) - Z,pj] dx = P,-(f, e), (8.6) 

JR JR 
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hence, by the resuhs of the previous section, it only remains to compute the right hand 
sides of the above expressions. 

We observe first that < C|(/)p hence we have 



N{(j))wj^xPj dx 



N{(j))ZjPj dx 



•7) 



cr' G {a, 3 /2a). To estimate the term involving C{(f>) we use integration by parts and 
the orthogonality condition (|4.6p to get 



^{(p)wj^xPj dx ^ I (f>[2wj^xxPj,x + Wj^xPj,xx] dx 

'^ziWj,xP])z + (f'{Wj.xPj)zz] dx 



+ / [/'(W) - f'iwj)](l)Wj,xPj dx 

^1 + 11 + III. 

To estimate the first term above we notice that in the support of pj^x and Pj^xx we have 



26 



25 



hence 



The second term is estimated similarly, using additionally = 0(a), so that 

\II\ < CQ,2-2a' + 15/16g-eoa|z|^ 

To estimate the last term we notice that 

|[/'(W)-/'K)K-|<C(a+l|5l + ||e||,,„), 

llGIlCG 

\III\<Ca'-''^'{a+\\S\\ + \\eU,^). 



To estimate the last term on the right hand side of the equation (|8.5p we observe 



that 



hence 



\w'{lj) - Wj,xPj\ < \w'{X,) - w'{x - f,)\p, + (1 - Pj)w'{lj) 
< Ca''\w'{x - fj)\p, + C(l - Pj)a^^^^^, 



E[w'iXj) - Wj^xPj] dx 



Analogous argument applies for the components of the right hand side of (j8.6p . Taking 
now cr < 2~^" we get that there exists a /i > 7/8 such that 



|Q,(f,e,^)| < Ca2+^ |P,(f,e,l)| < Ca'+f^. 



(8.9) 
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Now we will consider Lipschitz dependence of Qj and Pj on their parameters. To do 
this we can argue in a way similar to the way we proved Lemma 16.11 and Proposition 
16.11 This means wc need to introduce parameters g^. hj, f^, Sj. e^, fj and 6i and consider 
Lipschitz dependence with respect to them making use of Lemma 16.11 and Proposition 
16.11 In fact, when f is fixed then by a similar argument as above we get (using the 
familiar notation) 

liQf ^ - Qf ||e„« + ^P^'^ Pj'^We^c < Ca'^iM g'^\\e„a + ||h'/ - h^'||«„„) 

+ Ca'^ia\\{[-i;j\g,^ + \\e[-e'^\\eoc.) 

+ Ca'^dlei - e2\\e„» + a^WSi - (fall). 

(8.10) 

When all parameters except f are fixed then the argument is similar. Let us for example 
consider one typical term in the projection Qj: 

^''y^^Lf^l-^'^'y^^Lp^^d. < /;i0(^)-0(^)|kS.pg|dx 



+ ||fl -f2||oo||0<''||Moa,fJ 

<Ca2+Me-«o«|.|||f^_f^||^. 

Other terms can be handled in a similar way. In summary we will get then, taking into 
account the Lipschitz dependence on f only: 

IIQf - QflWa + ll^i'^ - P^,'^\\e„a < Ca^+^e-Oo-\^\\\{, - h\U (8.11) 

We will now consider the right hand sides of (|8.5p and ()8.6p . Let us go back to 
formulas (|7.4|) and (|7.6p . Using the same notation we will write now (|8.5p and (|8.6p in 
the form: 



- co/j' + 2VAra,o(5,(/j - /3,r,) sm{^ia,oz) + Cp(e"l/^"/-il - e-l^^-^^+^l) 

(8.12) 

(e;' + Xie,)do + hoj = He, + Pj, (8.13) 

where we have 

n^^. = n^^. (f " , e" , f ' , e' , f , ^) , fie, = fie, (f " , e" , f ' , e' , f , ^) . 

Functions Ily. and Ilg are smooth functions of their parameters and under the assump- 
tions jmil-jill) and ((3ll|) - ([3l4l) we have: 

||fi/JU„a + ||ne,||e„„ <Ca^ (8.14) 

The rest of this paper will be devoted to solving system (|8.12P " (|8.13p . 
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9 The Toda system 

In this section wc will determine the main order behavior of function f . Intuitively it 
should be given by solving the following system 

-co/j' + Cp(e-l-''^--^^-il-e-l-^^-^^+il) = 0, j = l,...,k. (9.1) 

Let us consider first the case k = 2. Remembering that in this case fs = oo and 
/o = — oo we can reduce (|9.ip to a single equation for u = fi — /2: 



(J 

u" + 2cpe"'' = 0, m(0)=: 01-02 Cp ^ , ^ „ , where C„ is defined in (frsl) . (9.2) 

and /j(0) = flj, see (|1.14p . For a > a family (parametrized by a) of explicit even 
solutions of (|9.2p is given by: 

1 , /-I +tanh^(Aaz/2)\ , , 

Ua(z) :=21og-+log2A-%-log( ^-^ ^J, a > 0, (9.3) 



where A = y ^I^-ai ■ Starting with this family it is natural to define 

fi{z) = -^u^{z), /2(z) ^ iM„(z). (9.4) 
Let us observe that in particular fi{z) + f2{z) = and 

\\f'^\U^O{a^), /3i=-^a, /2 - /i > 2 log - + log 2A2cp, (9.5) 

as needed. In this case we take 0o = and 9i = ^. In the sequel we will denote 

fo(z) = (/l(z),/2(^)). 

We will assume now fc > 2 since the case k = 2 has just been treated above. It is 
convenient to consider our problem in a slightly more general framework then that of 
the system (|9.ip . Thus for given functions qj{t),pj{t), j = 1, . . . ,k such that 

k k 



we define the Hamiltonian 



fe 2 



i=i j=i 



We consider the following Toda system 

dt 



dpj_ ^ _dH_ (9 6) 

dt dqj 
%(0)=9oj, Pj(0) = j = l,...,fc. 



Multiple-end solutions of autonomous elliptic problems 



Solutions to (|9.6p are of course even. Observe that also that the location of their center 
of mass remains fixed. Thus to mode out translations we will assume that 

k 

^go,=0. (9.7) 

We will now give a more precise of these solutions and in particular their asymptotic 
behavior as ^ — > ±00. To this end wc will often make use of classical results of Konstant 
PP] and in particular we will use the explicit formula for the solutions of (|9.6p (see 
formula (7.7.10) in gD]). 

We will first introduce some notation. Given numbers wi, . . . ,Wk G K such that 

k 

Wj = 0, and Wj > wj+i , j = 1, . . . , /c (9-8) 

we define the matrix 

Wo = diag(wi, . . . ,Wk). 
Furthermore, given numbers gi, . . . , gk G K such that 
k 

Y[gj = l, and5j>0, j = l,...,fc, (9.9) 

we define the matrix 

go = diag(gi, . . . ,5Ar). 

The matrices wq and go can be parameterized by introducing the following two sets of 
parameters 

= "'j - Wj+i, dj log.gj+1 - loggj, j = l,...,k. (9.10) 
Furthermore, we define fimctions <I>j(go, wq; i), t G M, j = 0, . . . , fc, by 

$0 = *fe ^ 1 

$j(go,wo;i)= (9.11) 
^ r,i...,^(wo)5ji . . .g^^. exp[-t(u),, + hu;^^.)], 

l<ii<-<ij<A; 

where r^^ . i^. (wq) arc rational functions of the entries of the matrix wq. It is proven in 
PO] that all solutions of (|9.6p are of the form 

gj(t) =log$j_i(go,wo;t)-log$j(go,wo;0, j = l,---,fc, (9.12) 

Namely, given initial conditions in (j9.6p there exist matrices wo and go satisfying (j9.8p - 
dnSl). According to Theorem 7.7.2 of [lOj, it holds 

(7j(+oo) = Wk+i-], 9j(-oo) = w-j, j = 1, k. (9.13) 
Wc introduce variables 

(9-14) 
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In terms of u = {ui, . . . , Uk-i) system (|9.6p becomes 

u"-7l/e" = 0, 
Uj{0) ^ qoj - qoj+i, Uj(0) = 0, j = l,...,fc-l, 

where 

/ 2 -1 O--- \ 
-1 2 -I--- 

7\/ = 

••• 2 -1 
V ••• -1 2 / 

As a consequence of (|9.1ip all solutions to (|9.15[) are given by 



(9.15) 



Uo{t) = qj{t) - qj+i{t) = 

Our first goal is to prove the following: 
Lemma 9.1 Let wq he such that 



-2 log $j (go , Wo ; < ) + log $j _ 1 (go , Wo ; i ) 
+ log$j+i(go,wo;t). (9.16) 



min [wj — Wj+i) — d > Q. 

j = l,...,k-l 



(9.17) 



Then there holds 

-Ck-jt - dk-j + r+(c) + 0(e-''l*l), as t +oo, j = 1, . . . , fc - 1, 



,it) = 



c-jt + dj +T^{c) + 0(e-''l*l), as t -oo, j = 1, . . . , A: - 1, 



(9.18) 

where T^ic) are smooth functions of the vector c = (ci, . . . , Ck-i)- 

Proof. This Lemma has been proved in [l^. We include a proof here for completeness. 
Let (jfj, J = 1, . . . , /c be a solution of the system ()9.6p depending on the (matrix valued) 
parameters Wq, go and defined in (|9.12p . We need to study the asymptotic behavior of 
4>j(wo, go; t) as i — > ±oo with the entries of Wq satisfying (|9.17p and still undetermined 
go- 

By (PH]) and we get that as i ^ -oo 

$j = (-l)^'('="^'Vi...j-(wo).9i . . ..gje-("'i+-+"'^)*(l + 0(e -(^^■-"'^-i)*) , 

hence 



$2 



^ .gj+i'^i...(j-i)(wo)ri...(,-+i)(wo)e ^ o(e-''l*l)) 



5j'^L.j(wo) 

It follows that as i ^ — oo 



(9.19) 



Uj{t) = log 



$1 



-Cjt + log 



= Cjt + dj+T^{c) + 0{e~''^'^) 



gj+i^i...(j-i)(wo)ri...(-,-+i)(wo) 
5j'^L.,(wo) 



+ 0(e-''l*l) 



(9.20) 
(9.21) 
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r?. ,(wo) 



where 

Similarly, as i — > +cx3 we get 



Uj (t) = log 



$2 



= -Cfc_,t-4-j+T+ + 0(e-''l*l), (9.22) 

where 



T+(c) = log 



?-fc+2-j...fc(wo)rfc_j-,,.fc(wo) 
'i+i-,...fe(wo) 



This ends the proof, q 



To find a family of solutions parameterized by a starting from a solution of (|9.6p we 



use functions uj and set 

Uja{z) — Uj{az) — 2 log — — log Cp. (9.23) 
Then functions fja{z) are obtained from the relations 

Uja{z) = fja{z) - fj+la{z), 

(9.24) 



Y.f,c{z)^0. 



Observe that as a consequence of Lemma l9. II we get that there exist Wj, gj, j = 1, . . . ,k 
such that (in;H])~(in3 holds, that 

min (wj — Wj^i) ~ ■& > 0, 

j=l,...,k 

and functions fj satisfy 

ll/,ll.^a = ll./^''"|^||loo<ca^ 

fj (z) ~ fj+i (z) > 2 log i + log Cp. 
In this case wc take 9o — and 0i — ^d. 



10 Bounded solvability of some equations in the line 

In this section will continue with preliminaries needed to solve (|8.12|) - (|8.13[) . First we 
will study the linearization of the system (|8.12p around the solution fo of the Toda 
system defined in (|9.4p and (|9.25p . We will always assume that a > is small and 
^0 > has value defined in the previous section. 
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Again, we will consider the case k = 2 first. Let Pj{z), j ~ 1,2 be given even and 
continuous functions. The linearized Toda system takes form 

ipi + Cpe"''°(v3i - 1P2) = Pi{z), 

ip'^-Cpe-''-{ipi-ip2)=P2{z), (10.1) 
^,(0) = x,, ^'(0) = 0, J = 1,2, 

which can be reduced to a single ODE for h = ipi — (p2 with p(z) = Pi{z) — P2{z): 

£„(/i) :=/i" + 2cpe-"°/i = p(z), z e M, 

h{0) = ho, h'{Q) = 0. ^ ' ' 

Notice that we do not impose the initial condition in (jlO.ip and (|10.2[) . In fact we will 
require that the solution to (jlO.ip and (jl0.2p is bounded and even and this determines 
uniquely the initial data. 

We will assume that function p{z) is even and satisfies 

Iblleoa < C'a2+'^, some/i>0, (10.3) 

and look for a solution to (|10.2p in the space of even functions such that 

\\h"\\0„a + a\\h'\\g„c. + a'^Moo < 00. (10.4) 

We will denote the space of such functions by X. Observe that evenness of solutions to 
(|10.2p is guaranteed by the initial conditions and the evenness of the right hand side. 

To solve problem (|10.2p we will construct a suitably bounded left inverse for the 
linearized operator By tjjja we will denote the elements of the fundamental set of 
£„. They are known explicitly: 

-01c(z) = W^(z), tp2a{z) = Zu'^{z) + 2, 

and their Wronskian is W{ipia, i^2a) = oP' ■ Then there exists a unique even and bounded 
solution solution to (|10.2p . It is given by: 

1 r 1 r 

h{z) = ^1pia{z) i'2a{s)p{s) ds + —^2a{z) ^/jia (s)p(s) ds 



a Jq a 

-^i'2a(z) 1 u'^{s)p{s)ds 

a Jo 

'^'aiz)u'^{s){z - S)p{s) ds + — / {u'^{s) - U'^{z))p{s) ds 



(10.5) 



1 

TV'2a(z) / u'^{s)p{s)ds. 

Directly examining this formula we get that 

\\h"\W» < C\\p\\g„,„ \hi^)\ < Ca-^\\p\\g,^ < Ca^, 

and using h'{Q) = we obtain the following estimate for the inverse of Ca'. 

\\h"\\eoo. + a\\h'\\g,^ + a^\\h\\^ < C\\pU,^. (10.6) 

These considerations provide us with a framework needed to solve (|10.2p . We have 
the following result. 
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Lemma 10.1 There exists a constants /i > 0, C > independent of a such that if p is 
an even function with 

iipIIsoo <c«'+^ (10.7) 

then Problem \10.2\ possesses an even solution h of the form HIO. 5\ such that the estimate 
ilO.6]) holds. Denoting this solution for a given p by TZa[p] we have additionally the 
following estimate 

\\T^a[Pl\" — 'na[P2\"\\eoa + aH'^abl]' - 'Tla[P2\' Wdaa + Oi \\Tla\pi\ — TZa [P'l] Woo 
< C\\pi -P2 Weaa, 

for all pi , p2 satisfying \10.7[ . 

The proof of this result is left to the reader. 

Now we will consider the general case k > 2. We are lead to the following linear 
system 



/ 2e 





V 



0... 



-e' 



\ 





XT 



(10.8) 



$= (^l,...,(/)fe_l), P= {pi,...,Pk-l), 

where p is an even function such that 



(10.9) 



We will analyze the solvability of this problem in the space of even functions </> such 
that 



||</'"||eoQ + Ol\\4>\\0oa + a^MWao < OO. 

Thus in addition to pO.Sp we will require that 

0(0) =.f, 0'(O) = O. 



(10.10) 



(10.11) 



We first observe that 



9j- 



dUr, 



a < 



1, j = m + 1, z OO, 

— 1, j = m, z ^ OO 

1, j = fc + 2 — m, z ^ — cx), 

— 1, j = k + I — m, z — OO, 



0, 



otherwise. 



Hence by a transformation we can find a set of linearly independent solutions to the 
homogenous version of (jlO.Sp 



Z OO, 
z —* — OO. 
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Similarly, considering derivatives of with respect to Wj we can find solutions of (jlO.Sp . 
4'2j{z), j = 1, . . . , fc — 1 such that 

V'2ja(t) =aej|z| + 0(l). 

The functions {{^ij{z), 'i{}[j{z)), (i/'2j(^), i^'2j{z))} form a fundamental set for the system 
(|10.8|) . whose behavior as 2; ^ ±00 is analogous to that of the functions {ilJi{z)ip'i{z)), 
(V'2(-z), V'2(^))i respectively. Let us denote the fundamental matrix of the system (|10.8p 
described above by ^'q(z) and the right hand side of the transformed system by q. We 
observe that as 2 — > 00, matrices 'i>~^ are block matrices of the form 

f al + oil) az/ + 0(l) \ ,_faI + oa) -azi + 0{1 



a/ + 0(1) y ' " \ aI + o{l) 

(10.12) 

Let us denote these blocks by ^mna, ^mna; rn,n — 1,2. respectively. Then, from 
variation of parameters formula we get that the solution of our problem has form 

{${z)J'{z))^^o,{z)- / ^-\s)-{G,q{s)Y ds^^,2o.{z) ^22Mq{s)ds. 
Jo Jo 

(10.13) 

Using this we can directly estimate 

form which it follows that (using the notation of Lemma 110.11 and going back to the 
original variables) 

mP\"\\eo^ + c^mPl'lKa + «'||7^[p1IU < C\\p\\e„a, 

(10.14) 

||7^[pl]" - n[p2]"\\eoa + a\\n[pi]' - n[p2]'\\eoa + - 7^[p2]|U < c\\pi -P2!le„a. 

(10.15) 

A second problem, important for our purposes is that of finding an even solution of 
the equation 

e" + Aie = g(z), z e M, (10.16) 
where, again, q is an even function with \\q\\a < +00. This time we want e to satisfy 

||e"||0oa + lle'llffoa + IWWdoa < +00. 

We need to assume the following solvability condition on q for this to be the case: 



q{z)cos{y/Xiz)dz = 0. (10.17) 

Under this assumption the solution turns out to be unique. Denoting the solution of 
(|10.16p by S[q] we get explicitly 



1 , — f°° — 1 , — f°° 

S[q] = sin(vAiz) / q{t) cos{\/ Xit) dt ^cos(vAiz) / q{t) sin{\/ Xit) dt . 

vAi Jz vAi Jz 

(10.18) 
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Clearly this operator is bounded in the sense that e = S[q] satisfies 

\\e\\9oc. + \\e'\\eoa + \\e"U,a < Ca-^\\q\\g„c.. (10.19) 

On the other hand, a better estimate is available in case that we know, in addition that 
Ik'UffoQ < +00. Indeed, in this case integrating by parts in the formula (jl0.18|) we get 

l|e||eoa + lle'lkoa + ||e"||e„a < C( ||q||(,„a + a-^\\q'\\9„c. )■ (10.20) 
Finally it is clear that the operator S[q] is Lipschitz in the norms used above. 

11 Solving the reduced system for (f,e,(5) 

We will now go back to (|8.12p - (|8.13|) . We will solve this system using a fixed point 
argument around and approximate solution f = fo, e = and 6 = 0, where fo is a 
solution to Toda system will be determined shortly. We let fo to be an even solution to 
(|9.25ll satisfying the initial data 

k 

/o,(0)=x„ 5]x,=0, 4(0) = 0, j = l,...,k. 

J=l 

We will set 

/3, =/o,(±oo). (11.1) 

In the sequel we will often use the fact that fg is a smooth function which in addition 
satisfies 

a-i|lClleoa + KWeoc + aWQeoc < Coa\ (11.2) 
With fo fixed now we will will assume that the solution of our problem is of the form 

f = fo+fi, 

where 

|!f;'||9„„ + a\\i[\\g,^ + a^lMoo < Ca^+\ (11.3) 

and for e and 5 we will assume that 

||e"||eoa + lle'llega + l|e||e„Q, < Ca^, 

(11.4) 

\5\ < Ca^+'', some z/ > 0. 
Let v he a fixed small number to be specified later and let 

||(fi,e,^lb = a-2(||friUoo + a||f{||^*„a + a'j|fi||oo) 1 



3^= |(fi,e,<5) 
Let also 



a-2+''(||e"||eoa + ||e'||eoa + ||e||eoa) + a"^ < cx) ^ 



BR = {if,,e,d) I ||(fi,e,<5)b <i?}. 
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We want to use Banach fixed point theorem to solve (|8.12p - (|8.13p in Bfja'^ , some i' > 0. 
It is convenient to linearize this system around the approximate solution first. We will 
denote fj = foj + fij and: 

L,(fi) = -ci/;; + Cp(e-(/«-^-/"^)(/y - - e-(^«^-^«^+^)(/i,+i - fij), 

- 2VAiajorfj(/j - /3jf?) sm{\/xlajQz), 
Kj(e) = [e] + Xiej)do, 

Soj = [^Oj,ii(fo) / Z"{s)Z{s)ds - XiAo,.22{io) [ Z\s) ds]cos{^ia,oz), 

hojifo + f,d)^ -/iOj(fo + f, rf) + djBoj. 

Let f , e and d e R'^ satisfying (|11.3p - (|11.4p be fixed. We will consider the linearized 
system written in the form: 

Lj(fi) Mj{f, d) + riy. (fo + f, e, d) + Qj(f , e, cT), (11.5) 
Kj(e) + JjBoj = S(fo + f , rf) + fie, (fo + f, e, d) + P,(f , e, d). (11.6) 

Our first goal is to show that given f , e and d one can chose 5 so that the solvability 
condition for the equation (|11.6p is satisfied. Thus we need to compute projections of 
(jll.6p onto cos(v^z) and show that adjusting 5 we can accomplish that for j = \, . . . ,k: 



Boj cos{^/ Xiz) dz ^0. (11.7) 



We first notice that given function q{z) such that 

""^Ik'lleoQ + IkllsoQ < oo, 

we have 



cos{^/Xl z) cos{^/XiajQz)q{z) dz ~ — / q{z) dz -\ — / coa{2 \/ Xiz)q{z) dz 



2./n ^ ' 2_ 



+ I zcos{y' Xiz)sm{y' Xiz)0{f3jr] )q{z)dz 
q{z) dz + O(a-l||q'||eoa) + 0{\\q\\g,,c.)- 



(11.8) 



The first term above which is of order 0(a~^ ||g||e(,Q,) is expected to be the leading order 
term. Thus to find asymptotic values of the integrals involved in (jll.7p we need to 
consider 



/o = / ^oj,ii(fo)rf2 / Z" {s)Z[s)ds, 

JQ JR 

/•oo /• 

//o = -Ai / ^oj,22(fo)d2 / Z^{s)ds. 

Jo JM 
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Let us recall that 

^Oj,ii(fo) = a,\((4f 

We will first prove the following 

Lemma 11.1 Let fo be the solution of the Toda system described above. For any B > 
there exists a cut- off function rj{t) =0, t < Ti, rj{t) ~ Q, t > T2 such that 

\h+Ih\>Ba. (11.9) 

Proof. Let us denote 61 = -j^{Z' f, 62 = -Xij^ Z^. Since a]^ = 1 + 0{f3]ri^), 
therefore 

j>00 

lo^bj ((4.)'-/3,V)rf^ + 0(a2) 

^0 

poo pOO 

= ^1 / ((4)' - /5f ) dz + biP^a-^ / (1 - v\t)) dt + 0{a^) 

JQ Jo 

= loi + Iq2 + Oia"^). 
Likewise we get 

IIo = b, r [((.4)' - P]) + WiViPi - /,')] dz 
Jo 

poo poo 

+ b2f3, / [(1 - r]^(t)){(3,a-^ - 1)] dt + 2(3ja-^ / ?7W('/W - 1) dt 
Jo Jo 

= //0l +//02 +//03- 

We notice that there exists a constant Co > 0, independent of the cut off function rj 
such that 

|/oi| + |//oi| <Coa. (11.10) 

To estimate the rest of the integrals we will consider 

b-s^biP^-a^ + (32l3ja-\ 

Notice that 63 = 0(1) as a ^ 0. Let us assume first that 63 7^ 0. Then we can take 
T2 = Ti -\- R (see the definition of 77) with R fixed. Increasing Ti if necessary we can 
make then |/o2 + as large as we wish while I//03I will remain bounded so that 

|/o + IIo\> \Io2 + //02I - I//03I - |/oi| - I//01I + 0{a^) > Ba. 

If on the other hand 63 = then by taking T2 — Ti = R large we can make //qs as large 
as we want and then 

|/o +//0I > I//03I - |/oi| - I//01I + 0(^2) > Ba. 

The proof of the Lemma is complete, q 

We are ready now to set up the fixed point argument. We will begin by outlining 
the steps that are needed to complete the proof. 
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1. Given (f,e, d) e 3^ we will chose S so that the solvability condition for pi.6p is 
satisfied. 

2. Using Step 1 we define an operator : Bjia" ^ 3^ by 

JF(f,e, d) — (fi,e, (5), where (fi,e, (5) is a solution of (|11.5|) - (|11.6p . 

3. We will show that there exist i/ > and R > such that 

WHleJ^yKRa-', (11.11) 

given that 

\\{leJ)\\y<Ra''. (11.12) 

4. We will show that the operator T restricted to -B_rq>^ is a contraction. 

Once steps (l)-(4) are executed will conclude by applying Banach fixed point theorem. 
Step 1 is an easy consequence of Lemma 111.11 and give the following formula 



S, I Bo, cosiVXiz) dz = / [hoj (fo + i,d) + n^, (fo + f , e, d) + (f , e, d)] 

X cos{\/\iz) dz. 



(11.13) 

We will need a precise estimate for S and this is the computation we will carry on now. 
We will make use of the following: 

Lemma 11.2 Let q : M_|_ —>■ M. be a function such that \\q'\\eoa < oo. Then for s e N 
we have 



q{z)z'^ cos( v \iz) dz 
q{z)z'' sm{\/Xiz) dz 



< 



c 



Yll^'llena, 



s > 0, 
s > 1. 



(11.14) 



Similar estimates hold when \\q Woga < oo, and and q'{0) ~ 0, with Wq'Wega above 
replaced by \\q"\\eoa- 

Proof. We have 

q{z) cos(-\/A]"z) dz = 



q'{z) sm{\/ Xiz) dz 



— I g"(z) cos(vAiz) dz, using (/'(O) = 0, 
^1 Jo 



hence we get 

q{z) cos{\/~\iz) dz 



< -Ik'lleoa, 
a 



q{z) cos( v \iz) dz 



<-h"\\eoo 
a 
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if ||f/'"||6ioa < oo- The first inequality in pi.l4p follows now by induetion. Similarly we 
get 

2:g(z) sin(\/A]"z) c?z = / cos{\/Xi z) dz — — / q' (z) sm{\/Xiz) dz 

vAi Jo ^1 Jo 



/ 2 1 \ f°° 

y-^JJ^ - —J J zq" {z) cos{\/\'iz) dz 



9 







2 . q (z) cos(\/Aiz) dz, 
'^i Jo 

using g'(0) = 0. 

From these identities, using the above and induction, one can show the second estimate 

in mniD. □ 

With some abuse of notation we will write 

|l(f,0,0)b = !|fb, ||(0,e,0)b = |leb, etc. 

Will will now compute the integrals on the right hand side of (|11.13p . Wc begin 
with 



[/iOj(fo + f)] cos(vA7z) dz 

w"{s)Z[s)ds I AQj iiajQCOs{\/\iz)dz 
i Jo 

— 2dj / (Z'(s))^(is / AQj^i2ajQy/}^sm{y/J^ajoz) coa{y/\iz) dz 
Jm Jo 

+ 4 / + /j - f],7]\\eoa + a2)e-''°"l^l cos{y%z) dz 

Jo 

I w"{s)Z{s)ds)lj -2dJ I {Z'{s)f ds\llj+dpllj. 
Since fo is a smooth function, therefore from (|11.2p we get using /jo(0) = and Lemma 



Ij = /_ %o ((4- +f-f- Ppf) cos( vOTz) dz 
%o((4)'-/3'^7'))cos(v^z)dz 
+ / a,o(2/;o/i + (/;)') 003(7X7) dz, 







hence 



Similarly we get 



\I,\<C{a'+a^f\\y). 

\II,\<Ca\\d\\y{a^ +a\\{\\y), 
\III,\<Ca^dfyil + a\\{\\y). 
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Some tedious but standard calculations yield: 



Ilej (fo + f , e, d) cos( V Aiz) dz 



< C[a^ + a\\d\\y{a^ + a\\i\\y + a^-^My) 



Using (|8.8p and Z{s) « e "2 ^ I'*! , p > 2 we can estimate 



Pj (f , e, d) cos( v Aiz) dz 



Since a' < 2 the above estimates and Lemma [11.11 yield whenever ly < j'. 

\\S\\y<^[ai+R'a''], 



(11.15) 



as long as (|ll.lip is satisfied. 

Step 2 follows directly from the results of Section [TOl We will now show that the 
claim of Step 3. From (jlO.lSp we have 

e, = -5jS[Boj] + 5[S-(fo + f , d)] + S[Il,^ (fo + f , e, d)] + 5'[P,(f , e, d)], (11.16) 
where S is the operator defined in (|10.18p . Wc have, using estimates p0.19p - p0.20p : 

\\S,S[BoMeoc.<Ca^S\\y. 

Likewise 

\\S[hoj{io + I rf)]||eoa < i?oa' + Calmly + Ca^MyfiWy. 

where Rq > is a constant depending on fo only. Using the formulas for He- and Pj 
we get 

\\S[Ile^{io + le,d)]\\e,^ < Ca-^[a^--'\\e\\y{a\\d\\y + a^)], 
\\S[P,{le,d)]\\g,^<Ca^+i, 

These estimates can be summarized as follows: 

\\e\\y < Roa" + Ca''{a^ +0"). (11.17) 

We will now estimate the right hand side of pi.5p . 

||A4j(f,dll|eo« <^0«'l|rf||j^ + C«2(||fi|2^ + ||dj|^||f||^). 

We also get in Bjia-^ : 

WUf^ifo + le, d)\\g,^ + \\Q, ile, d)\\g,^ <Ca\a-s +a^''). 
The last two estimates yield 

llfib < Ro\\d\\y + C{a^ + a^n- (11-18) 
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Combining now (|11.15p - (|11.18p we get 

\\J^{le,d)\\y < Roa-' + ^^^^^-^{ai + a-') + C{ai + a^'') (11.19) 

< f (11.20) 

provided that ^ < ly < j, R in the definition of B^a" is taken sufficiently large, a is 
taken small and B in the Lemma [11.11 is taken large. This shows the claim of Step 3. 

It now remain to show the claim of Step 4. This follows from a rather straightforward 
application of the theory of Lipschitz dependence of various terms involved in (jll.Sp - 
(|11.6p . We have for instance 

||^^(^)-5^(^)b<§||(f(^\e(^\<?^))-(f(^\e(2\dt^))b, (11.21) 

and 

We leave the details to the reader. From these estimates, taking B larger if necessary 
we get that is a Lipschitz contraction as claimed. This shows that T has a fixed 
point in Bna" . The proof of the theorem is complete, q 



Acknowledgments: This work has been partly supported by Chilean research 
grants Fondecyt 1070389, 1050311, FONDAP, Nucleus Millennium grant P04-069-F, 
an Ecos-Conicyt contract and an Earmarked Grant from RGC of Hong Kong. The 
fourth author thanks Professor Fanhua Lin for a nice conversation on the classification 
of solutions to (|l.ip . which motivated this research. 

References 

[1] Ambrosetti, A., Badiale, M., Cingolani, S., Semiclassical states of nonlinear 
Schrodinger equations. Arch. Rational Mech. Anal. 140 (1997), 285-300. 

[2] Byeon, J., Wang, Z.-Q. Standing waves with a critical frequency for nonlinear 
Schrodinger equations. Arch. Rat. Mech. Anal. 65(2002), 295-316. 

[3] H. Berestycki, P.-L. Lions, Nonlinear scalar field equations. I. Existence of a 
ground state. Arch. Rational Mech. Anal. 82 (1983), no. 4, 313-345. 

[4] J. Busca, P. Felmer, Qualitative properties of some bounded positive solutions 
to scalar field equations. CaZc. Var. Partial Differential Equations 13 (2001), 
no. 2, 191-211. 

[5] V. Coti Zelati, P.H. Rabinowitz, Homoclinic type solutions for a semilinear 
elliptic PDE on R'\ Comm. Pure Appl. Math. 45 (1992), no. 10, 1217-1269. 

[6] M. Crandall, P. Rabinowitz, Bifurcation from simple eigenvalues, J. Funct. 
Anal. 8 (1971), 321-340. 

[7] E.N. Dancer, New solutions of equations on R". Ann. Scuola Norm. Sup. Pisa 
CI. Sci. (4) 30 no. 3-4, 535-563 (2002). 



Multiple-end solutions of autonomous elliptic problems 



[8] C. Delaunay, Sur la surface de revolution dont la courbure moyenne est con- 
stante, Jour, de Mathematiques 6 (1841), 309-320. 

[9] M. del Pino, P. Felmer, Local mountain passes for semilinear elliptic problems 
in unbounded domains, em Calc. Var. Partial Differential Equations 4 (1996), 
121-137. 

[10] M. del Pino, P. Felmer, Semi-classcal states for nonlinear Schrodinger equa- 
tions, J. Fund. Anal. 149(1997), 245-265. 

[11] M. del Pino, M. Kowalczyk, J. Wei, Concentration on curves for nonlinear 
Schrodinger equations, em Comm. Pure Appl. Math. 70(2007), 113-146. 

[12] M. del Pino, M. Kowalczyk, J. Wei, The Toda system and clustering interfaces 
in the AUcn-Cahn equation. Arch. Rat. Mech. Anal., to appear. 

[13] Floor, A., Weinstcin, A., Nonspreading wave packets for the cubic Schrodinger 
equation with a bounded potential, J. Fund. Anal. 69 (1986), 397-408. 

[14] B. Gidas, W. M. Ni and L. Nircnberg, Symmetry of positive solutions of 
nonlinear elliptic equations in , Adv. Math. Suppl. Stud. 7A (1981) 369- 
402. 

[15] A. Giercr and H. Meinhardt, A theory of biological pattern formation, Kyber- 
netik (Berlin) 12 (1972) 30-39. 

[16] P. Gray and S.K. Scott, Autocatalytic reactions in the isothermal, continuous 
stirred tank reactor: isolas and other forms of multistability. Chem. Eng. Sci. 
38 (1983), 29-43. 

[17] C. Gui, J. Wei, Multiple interior peak solutions for some singular perturbation 
problems, J. Differential Equations 158 (1999) 1-27. 

[18] Kang, X., Wei, J. On interacting bumps of semi-classical states of nonlinear 
Schrodinger equations.. Adv. Diff. Eqn. 5(7-9), 899-928 (2000). 

[19] N. Kapouleas, Complete constant mean curvature surfaces in Euclidean three- 
space, Ann. of Math. (2) 131 (1990) 239-330. 

[20] B. Konstant, The solution to a generalized Toda lattice and Representation 
Theory, Adv. Math. 34 (1979), 195-338. 

[21] M.K. Kwong, Uniqueness of positive solutions of Am — u + u^ = in i?". Arch. 
Rational Mech. Anal. 105 (1989), no. 3, 243-266. 

[22] Li, Y.Y. On a singularly perturbed elliptic equation, Adv. Diff. Eqn. 2(6), 
955-980 (1997). 

[23] F.H. Lin, W.-M. Ni, J. Wei, On the number of interior peak solutions for a 
singularly perturbed Neumann problem. Comm. Pure Appl. Math. 60 (2007), 
no. 2, 252-281. 

[24] A. Malchiodi, Some new entire solutions of semilinear elliptic equations on 
i?". Preprint 2007. 

[25] A. Malchiodi, M. Montenegro, Boundary concentration phenomena for a sin- 
gularly perturbed elliptic problem. Comm. Pure Appl. Math. 55 (2002), no. 
12, 1507-1568. 

[26] A. Malchiodi, M. Montenegro, Multidimensional boundary layers for a singu- 
larly perturbed Neumann problem.Dufce Math. J. 124 (2004), no. 1, 105-143. 



Multiple-end solutions of autonomous elliptic problems 



[27] R. Mazzeo and F. Pacard, Constant scalar curvature metrics with isolated 
singularities. Duke Math. J. 99 (1999), no. 3, 353-418. 

[28] R. Mazzeo and F. Pacard, Constant mean curvature surfaces with Dclaunay 
ends. Comm. Anal. Geom. 9 (2001), no. 1, 169-237. 

[29] R. Mazzeo, F. Pacard and D. Pollack, Connected sums of constant mean 
curvature surfaces in Euclidean 3 space. J. Reine Angew. Math. 536 (2001), 
115-165. 

[30] R. Mazzeo, and D. Pollack, Gluing and moduli for noncompact geometric prob- 
lems. 17-51, Synipos. Math., XXXVIII, Cambridge Univ. Press, Cambridge, 
1998. 

[31] R. Mazzeo, D. Pollack, and K. Uhlenbeck, Connected sum constructions for 
constant scalar curvature metrics, Topol. Methods Nonlinear Anal. 6 (1995), 
207-233. 

[32] W.-M. Ni, I. Takagi, On the shape of least-energy solutions to a semilinear 
Neumann problem. Gomm. Pure Appl. Math. AA (1991), no. 7, 819-851. 

[33] W.-M. Ni. I. Takagi, Locating the peaks of least-energy solutions to a semi- 
linear Neumann problem. Duke Math. J. 70 (1993), no. 2, 247-281. 

[34] W.-M. Ni. Diffusion, cross-diffusion, and their spike-layer steady states. No- 
tices Amer. Math. Soc. 45 (1998), no. 1, 9-18. 

[35] Oh, Y.-G. Stability of semiclassical bound states of nonlinear Schroedinger 
equations with potentials, Commun. Math. Phys. 121 (1989), 11-33. 

[36] P. Quittner, P. Polacik and Ph. Souplet, Singularity and decay estimates in 
superlinear problems via Liouville-type theorems, I: Elliptic equations and 
systems, Duke Math. J. 139 (2007), 555-579 

[37] R. Schoen, The existence of weak solutions with prescribed singular behavior 
for a conformally invariant scalar equation. Comm. Pure Appl. Math. 41 (1988) 
317-392. 

[38] Strauss, Walter A. Existence of solitary waves in higher dimensions. Gomm. 
Math. Phys. 55 (1977), no. 2, 149-162 

[39] A. M. Turing, The chemical basis of morphogenesis, Phil. Trans. Roy. Soc. 
Lond. B 237 (1952) 37-72. 



